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ABSTRACT 

This presentation sets forth a method for computing the 
three-dimensional transonic flow around the blades of a compressor or 
of a propel leTi The method is based on the use of the velocity 
potential, on the hypothesis that the flow Is Invlscid, Irrotatlonal 
and Isentroplc. 

The equation of the potential is solved in a transformed space 
such that the surface of the blade is mapped Into a plane where the 
periodicity is Implicit. This equation is In a nonconservative form 
and is solved with the help of a finite difference method using 
artificial viscosity and artificial time. 

A computer code is provided and some sample results are given in 
order to demonstrate the influence of three-dimensional effects and the 


blade's rotation 


I. INTRODUCTION 


Many scientists and engineers continue to study the ways and means 
of making better use of the energy resources at our disposal, even 
though the energy crisis is perhaps no longer considered so severe. 
Tlirbines and compressors, which both create and consume energy, have in 
recent years been the subject of many theoretical and experimental 
studies aimed at improving their design and efficiency. Although some 
major Improvements have been introduced lately, the study of the 
transonic flow across a single stage of a turbine, or of a compressor, 
still remains highly complex. For this reason most of the theoretical 
work done so far has focused on the two-dimensional cascade problem 
[13-17] or on the mean flow problem [22]. 

Let us first discuss the physical background. A basic compressor 
consists of a succession of rotors and stators. These rotors and 
stators are situated on the hub and surrounded by the cowling; they are 
composed of a certain number of blades distributed around the hub, the 
shape of the blades depending on their use. If the hub is cut along a 
generator line and transformed into a plane, a so-called ’’cascade" of 
these blades is obtained. A propeller is a compressor without cowling. 

The present study deals with the problem of transonic flow around 
compressor or propeller blades. From the mathematical point of view, 
this leads us to a system of partial differential equations of mixed 
type, in which the unknowns are the geometric and physical 

characteristics of the compressor. For a three-dimensional analysis 
these equations are too complex to be Integrated without some 

simplifications. First of all, we suppose that the fluid we are 
concerned with is a poly tropic and nonvlscous gas and that a velocity 


potential exists. This necessitates a further hypothesis, namely, that 
the variation of the entropy is small so that the entropy itself 
remains essentially constant. Thus our system of partial differential 
equations becomes equivalent to a single equation, the potential 
equation. In both the steady and the time-dependent cases. 

The potential equation has already been solved numerically in the 
three-dimensional case for oblique and swept-wings [2,14,15]. However, 
because of the periodicity of the compressor problem, the square-root 
transformation used in those works is not practical here. We therefore 
propose a new transformation, which maps the surface of the blade into 
a plane and Includes periodicity implicitly. 

The scheme we use Is similar to the one used by Jameson and 
Caughey [15] in the development of the swept-wlng code known as FL022. 
We solve the finlte-di f fere nee approximation to the potential equation 
by row relaxation. A typical run consists of 50 iterations on a 48^6x4 
grid, followed by 1 00 iterations on a 96xl3<8 refined grid. This takes 
15 minutes on the GDC 6600 and 3 minutes on the STAR. The simplicity 
of bur grid generation, together with the other hypotheses given above, 
limits our choice of blade geometries. Nevertheless, this method 
enables us to study how the speed of rotation influences the relative 
flow around the blades, and we have compared our three-dimensional 
results with data for two-dimensional cascades. 

In Section 2 we shall derive the equations of motion In physical 
space and, in Section 3, we shall consider the potential equation, 
which is obtained after several changes of variable. The numerical 
scheme Is presented in Section 4. The numerical results obtained are 
shown in Section 5 together with some Calcomp plots, while Section 6 
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provldes a maiual on ttie um of our computer uode^ Sections 7 and 8 
contain a bibliography and a listing of the code. 

This study has been supported by NASA under Grant No. 
NGR-33-016-201 and by tlie U. S. Department of Energy under Contract 
DE-AC-02-76EW)3077 and I take this opportunity of expressing my sincere 
gratitude to Prof. P. Garabedian for his invaluable advice and to Dr. 
F. Bauer for her constant encouragement. I am also indebted to my 
entire family for their moral support. 
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II. THE PARTIAL DtpreRENTIAL EQUATIONS OF MOTION 
This section sets forth the various equations used in our 
computational method and in the resulting code. 

The general equations of fluid dynamics for an inviscid gas are 
well known [6). These are: 

(a) the equation of conservation of mass, 


(1) P(. + pUj^ + pVy + pW^ + PjjU + PyV + PjjW » 0, 

where p Is the density, (u,v,w) are the velocity components, and 
(x, y,z) are the coordinates In physical space; 

(b) the equations of conservation of momentum, 

p(u^ + uUjj ■¥ vuy + wu^) + Pjj » 0 

(2) p(v^ + uVjj -f- vvy + wv^) Py » 0 

p(w^ -f UWjj + VWy + WWjj) + P^ =* 0 , 

where P is the pressure; and 

(c) the equation of conservation of energy, 

( 3 ) §- 0 . 


where S is the entropy. It will be recalled that 


d 

dt 



^ + v ~ + w ^ 
3x 3y 3z 


is the naterial time derivative. We are neglecting here such external 
actions as gravity. 

For the following calculations we shall assume that the entropy is 
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constant throughout the fluid; the last equation will therefore not be 
used. But then, In order to represent weak shock waves which occur In 
transonic flows, we have to replace the Ranklne Hugoniot shock 
conditions. This Is done by permitting a jump In the horizontal 
component of momentum and by adding artificial viscosity terms to the 
partial differential equations. The approximation is adequate for Mach 
numbers close to 1. 

Let us suppose now t}»t the fluid is a polytropic gas. Therein: 
(4) P - A(S)pY . 

where the function A(S), on the basis of our hypothesis, becomes a 
constant, and where 




is the adiabatic exponent of the gas, l.e. the ratio of the specific 
heats at constant pressure and at constant volume. 

If c Is the local speed of sound In the gas, we have 


(5) 



so that equations (2) become 

p(Ug -f UU^ + VUy + WUg) + c^p^ * 0 , 

p(vj + uVjj + vvy + wVjj) + c^Py - 0 , 

p(Wj. H- UWjj + VWy + WW^) + c Pg ■ 0 . 


( 6 ) 
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Suppose now tliat the flow Is steady, Ite. that the flow Is independent 
of time, so tfiat 



Then from (6) and (7) we obtain 


(S) P(ll^ + Vy + Wjj) - (U^U^+ UVUy+ UWUjj 

2 2 

+ VUVj^+ V Vy+ yWY^"^ WuWj^+ WVWy+ W Wjj) * 0 s 

or, if we assume the existence of a velocity potential, (^ , and If we 
divide by p /c^, 

( 9 ) c (♦ JJJJ+ yy+ ♦ 2 2 ^ “ xx*^ y y^ ^ z 2 

+ 2uv(j»jjy + 2vw(py^-h 2wu)>2jj) - 0 , 

with u » (p^, V » ^y, w ■ (^ 2 * potential equation. It is 

hyperbolic for supersonic flow, 

u^ + w^ > , 

and elliptic for subsonic flow, 

u^ -f + w^ < c^ . 

Equations (4) and (5) show us that 



( 10 ) 


I d dp 

FT'-J"-* 
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Using tlie velocity potential, equations (6) become 


(^(♦J+,2+,2) +^)^. 0 . 

(11) (♦$ +*y +*^) +^)y - 0 . 

(7 (*2 . 0 ^ 


so that 

2 

(12) i + 4>y + 4)^) + » constant 


throughout the fluid. This Is the Bernoulli equation, which enables 
us, knowing the velocity potential, to compute the speed of sound. 

Suppose next that the flow is no longer steady with respect to the 
initial frame of reference, but that. If we consider a frame turning 
around the x-axis with a constant speed of rotation, o) , the flow is 
again Independent of time. Then, Instead of (7), we have 


3p j. 

■ - (t)Z + wv — . — 
3t 3y ^ 3z ’ 3t 


3u . 3u 
noz o)y ~ 

3y ■ 3z 


3v 3v. 3v 3w 

- OJZ — T 0)V — . — 

3t 3y ^ 3z * 3t 


3w . 3w 

- (uz ~ + uy — 
3y 3 z 


( 13 ) 
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By introducing the cylindrical coordinates (x,9,R), the equations (13) 
become 


i£- 

3t 


3p 3u ^ 3u 3v ^ _ 3v 3w 
3 ? ’ Ft " ” ^ 39 ^ 3 1 " “ “ 30 ’ 3 1 



These new equations show us that the flow dependence on t and 9 is 
characterized by a dependence on (9 - ut) alone* 

In this time-dependent case a combination of equations (1) and (6) 
gives us 

(14) P(U^ + Uy -f Wjj) - (ugu^+ UqVoU + UoWqUj, -f VqUqVx + v§Vy 

where components of the relative velocity defined by 

(15) Uq ■ u , vq ■ V - uz , wq - w + wy . 

If we again assime the existence of a velocity potential. It must now 
satisfy the equation, 

(16) (^^xx+ VQf^yy + w§^2 2‘*' ^UQVo^xy'*' 2 wquq^zx^ " ° 

This equation is similar to equation (9); it is a second-order 
nonlinear partial differential equation which is hyperbolic for 

u§ + v§ + w§ > c2 , 
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and elliptic for 



11 ^ 4 - 4 . rj2 y _2 

Uq T Vq T Wq \ C 


We still have the relation (10) for the speed of sound, so that, 
using the velocity potential and with the help of (13), we obtain, 
instead of (6) , 


(17) 




and this gives us the new Bernoulli equation, 

2 

(18) J (4'x‘^y'N>z) ■* wz4>y + d- — ■ » Constant 

along each line parallel to the x-axis. 

In the next section we shall transform these equations, by a 
change of coordinates. Into the system which is solved by our code. 
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III. GRID GENERATION 

We have obcalned the equations to be solved, but we have still to 
Impose the geometrical constraints due to the slip condition on the 
surface of the compressor blade and to the periodicity condition. In 
this section, we shall describe the geometrical space In which the 
equations will be solved. We shall therefore list all the mappings 
which are performed to transform the physical space onto the 
computational space. 

For the purpose of the periodicity condition, we begin by 
Introducing the angle 0 of the cylindrical coordinates (x,0 ,R) with 
respect to the x-axis. To accomplish this we use a conformal mapping 

(19) (x,y,z) -► (xo,0,Zq) 
defined by 

(z+ly) “ exp(ZQ+l0), x » Xq . 

The Jacobian natrlx associated with this transformation Is 

10 0 

(20) * 0 P Q , 

0 -q P 

with elements defined by 

g cos6 ^ sln9 N 

P - — g — , q - — g — . R ■ exp(ZQ) . 


We also require the Inverse matrix, namely. 
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1 0 0 

■ 0 z -y . 

0 y z 

If (t,j,k) is the orthonormal basis for physical space, then the basis 

> > 

vectors connected with the new coordinates are I, u, and v, where 

•► +• + 

(21) u - 2 j - yk , V « yj + zk . 

For the surface of the blade, the representation of the finite 
difference slip condition becomes greatly simplified and more accurate 
If the boundary surface lies on a coordinate plane. The idea, cf. 
[lA], is to transform the surface of the blade onto a plane which will 
constitute the lower boundary of a half -space. 

Unfortunately, if we apply the sq\iare-root transformation of 
reference [14], the periodic strip is transformed in such a way tJiat 
the periodicity condition is difficult to satisfy. We therefore prefer 
a transformation which would map a periodic strip conformally onto a 
half-space, so that the periodicity condition becomes implicit. 

The required transformation can be decomposed into two successive 
mappings. For the putpose of simplification we consider only plane 
sections orthogonal to the Z^-axls. The first mapping transforms a 
periodic strip (cf. Figure la) onto a slit plane (cf. Figure lb). 
The image of the two lines delimiting the strip is the negative real 
axis. We then apply a square-root transformation which maps this plane 
onto a half-plane (cf. Figure Ic). 



Figure la. The periodic strip 



Figure lb. The slit plane. 






To Implement these transformations we draw in each plane section a 
singular line (the branch cut) from Inside the blade, near the nose, 
out to downstream Infinity. This is actually a half-line whose origin 
has the coordinates (Xg,9g) . We next perform the mapping 

(22) (xq,9 .Zq) (X,Y,Z) , 
defined by 

2 

(Xq-x,) + 1(9-«8) - N"‘ Xog( 1 + , Zfl * ^ • 

where N and T are two constants. N Is more precisely the number of 
blades on the compressor. The coordinates Xg and 9g usually depend on 
Z. Their derivatives with respect to Z will be denoted x^ and 9^ , 
respectively . 

The Jacobian matrix determined by this transformation is 
a -b 0 

(23) J 2 - b a 0 , 

e f 1 

whose elements are defined by 

2 2 

a ■ Hxj^ , b ■ H9jj , e ■ - ax^ ” b9 2 t f “ “ "*■ ^ ®X * 


The Inverse matrix Is 
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^ X ^ 

•^2^ - “®X ^ 

3 7^ I 

and the basis vectors related to 
with 


OF 


» 


this transformation are now (A,B,C), 


+ •*■*•**• **■ 

(24) A “ Xjji + 6 j^u , 8 = - 0^1 + x^u , C“X2l+3;2U + v 

Let ns return to the potential equation (9), tt can he 
refornulatad as 

(2 5) - (7^ . V)2^ = 0 , 

where 


V = -Li + i_j4--Lk 

3x 3y 3z 


is a notation for the gradient • If we denote the final Jacobian matrix 


A » J j • J 2 


by 


and, if we use 


(aij) 


-4-3 .4 3 . 3 

1 “ 3X * 2 * 3Y * *^3 “ az 



0RlGir>3AU 1^* 
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for the derivatives with respect to the new coordinates, then the 
derivatives with respect to the Initial coordinates are represented by 

h " ' h' \ ' 

Thus the Laplactan, which is defined by 



can be rewritten as 

(26) v 2 - I + I a^ja^^d,^ , 

i*j,k ■ l.j.k 

where the coefficients a^j stand for the derivatives of the Jacobian 
matrix elements; 

®lk • 

Let us set 

(F^j) » F - A*^*A , = 7^X , D 2 = 7^Y , D 3 = V^Z . 

These notations allow us to simplify the equation (26) to arrive at 

(27) 72 - I Fj k ijlk + I Mk • 

j,k k 

For the second term of equation (25) let us denote the components 
of the velocity in the physical space by 
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3^ 


3(^ 


3^ 


. * 3-11 


Using the same computation as for the Laplaclani we obtain 


(28) (V<t>.V )2 - I + y Cj^djj , 

k,£ I 


■¥ •¥ * 


where ( 8 ^ 82 , 83 ), the components of the velocity In the basis (A,B,C), 
are given by 


® j * I ^i I I 


and where the coefficients Cj^ are defined by 

= (7^»7)h , C 2 = (7^»7)^ y, C 3 » (V(j)«7)2 z . 

After these necessary but somewhat tedious transformations the 
potential equation acquires the useful form 


(29) I (c^F. . - B.Bjd.d.* + [ (c^Dj ^ Ci) » 0 

IJ i 


This has the advantage of being relatively tractable, for an equation 
which Is, really, quite complicated. In the time-dependent case, we 
need only to replace the <^j's by the components of the relative 
velocity in the calculation of the coefficients of equation (29), 
Moreover, if a reduced potential Is defined as the difference between 
the true potential and the uniform flow potential corresponding to the 


inlet speed* we can use it to compute the derivatives dj^dj^ and d^^ in 
equation (29) without changing that equation. 

We next describe the grid used in our computational method. It is 
obtained with the help of a system of sheared coordinates. Itiese are 
defined by considering coordinates parallel to the transformed surface 
of the blade. If 

Y - S(X,Z) 

is the equation of this surface* the transformation in question is 
defined by setting 

I » X , J - Y - S(X,Z) * K - Z . 

The Jacobian matrix is 

I -Sx 0 

(30) J3 - 01 0 , 

0 -S2 1 

and the related basis vectors are 

(31) I - A + S^B , J - B , K « C + S28 . 

In order to carry out the computations on a finite domain, these 
last coordinates are stretched so that the final domain of computation 
becomes a cube the edges of which have length 1. After these 
transformations equation (29) remains of the same general type* but its 
coefficients have to be changed slightly. Complete details are 

specified in the listing of the code in Section VIII. 
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It remains to define the boundary conditions at the hub, at the 
cowling, arel on the blade surface. We Impose a slip condition for the 
flow on these boundaries; this leads to a Neumann problem for the 
velocity potential, since the normal derivatives at the boundaries must 
vanish. We shall express this condition In our system of coordinates. 
The huh and tlie cowling are defined by equations of the form 

7 , ■ Constant . 

At any boundary points the two tangential vectors are I and u, and v ^ 
i X u is a normal vector. The boundary condition at the hub and at the 
cowling may thus be expressed by the equation 

(32) ^ = y.^2 + ztl * 0 • 

3Zo 

On the other hand, the blade surface Is defined by the equation 
J * 0 . 

> > •> > 

The two tangent vectors I and K lead to the normal vector I x ’< whose 

coordinates In the basis (i,u,v) are 

- Sj^e — f h Sr; 

c =■ R-=(9x + Sj(Xx) . 9 - Sx9x-*x • ’f ‘ ff ^ • 

Hence the boundary condition on the blade surface Is given by 
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OS) a ^ + 0 |i +Y Yy)t2 ‘*' • 0 . 


In equations (32) and (33), (<^i,^2»'^3^ derivatives of 

the true potential, i.e* the components of the velocity. For the 
time-dependent case, they are replaced, in equation (33), by the 
components of the relative velocity. 

In the code the values of (<^pi^2*'^3^ computed only when 

necessary, and the only derivatives available are 



3^ 

n 


W 


3 ^ 

T? • 


These derivatives are related to the components of the velocity through 
the relation 

<^1 U 

(34) ^2 • ^ ^ • 

♦ 3 W , 

The boundary conditions (32) and (33) become 

(35) eU + fV + W - 0 
and 

(36) (oa + 0b + ye)U + (-ab + 0a + yf)V + yW - 0 . 

In three dimensional space an appropriate vortex sheet behind a 
blade niist be considered. The shape of the vortex sheet is modeled by 
our code in the following way. In each plane cross-section, the upper 
and lower surfaces of the blade are extended behind the trailing edge 
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by two lines parr<llel to a branch cut (Cf. Figures 1,4,7). These two 
lines represent the vortex sheet in that section, and If the blade is 
closed they are identical. For two points situated on opposite sides 
of the vortex sheet the pressures are the same, and the normal 
velocities are zero; only the tangential components of the velocity may 
be different. Since a shape is assumed for the vortex sheet, we 
require only continuity of the normal component of the velocity across 
the sheet. Conputationally this reduces to a condition like 

( 37 ) ® 

after a jump in ^ is removed. Moreover, the jump of the potential ^ 
across the vortex sheet Is supposed to be constant in each plane 
section. The Kutta-Joukowsky condition at the trailing edge determines 
this jump. This amounts to a linearized treatment of the vortex sheet. 
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IV. FINITB DIFFERENCE SCHEME 

the success of codes for the design and analysis of supercritical 
wings [1,2,3,15] shows how effective the computHtlona I fluid dynamics 
has become for transonic flow. The first step In this development was 
the design of shockless airfoils by the hodograph method [{]• Then the 
introduction of a retarded difference scheme [20] allowed the analysis 
of flow at off-deslgn conditions. This scheme incorporated artificial 
viscosity in order to capture shocks In the supersonic zone. It was 
then Improved to permit the analysis of three-dimensional wings [14] 
and their design [ll]. In our computational method we use the 
last-mentioned scheme, and this section explains how it Is incorporated 
Into our computer code. 

All the equations we have described In the previous section are 
represented In the cpraputer code by finite difference approximations. 
The flow conditions at each grid point are determined with the help of 
the reduced potential G. The first derivatives of this potential are 
calculated by using central differences. With these values we compute 
the approximate velocity. This allows us to determine, with the help 
of Bernoulli's equation (12), whether, at the point considered, the 
flow is subsonic or supersonic^ 

If the flow is subsonic the second derivatives used for the 
computation of the potential equation residual are approximated by 
central differences of the form 

. _ ^^i+l,j,k - + Gi-i,j,)c 

Gxx - —2 . 


( 38 ) 



and 


" 23 “ 


ORIGINAL FASl W- 
OF POOR QUAU’^ 


(39) G 


XY 


°l+l, j+l,k “ j-i,k " ^l-l,J+l,k 


For supersonic points, on the other hand, we have to introduce 
artificial viscosity in order to capture weak shocks. This Is 
accomplished by using a retarded difference scheme [13]. Thus we 
separate the equation (29) into two groups of terms: 

(40) (c2-q2) Ggg + c2(v2g - G^g) =» 0 , 


where q is the speed and s a coordinate in the flow direction. The 
first term represents the second derivatives in the flow direction; tlie 
derivatives in tl»e other directions form the second term. The second 
derivatives used in this second expression are computed by using 
equations (38) and (39). But for tlie second derivatives of the first 
expression we use retarded differe :e approximations of the following 
types: 


(AX) 2 


and 


(4 2) r. 


XY 


_ ^'l,.1,k ■ ^’l,j-l,k - ^l-l,j,k +/’n-l,j-l,k 


AX AY 


provided that the velocity luis positive components in the X direction. 
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and in the Y direction. Equations (41) and (42) are only first-order 
accurate and Introduce the truncation errors 


■ ^X^XXX • 


AXG^XY ®XYY 
5 


For the potential equation (40), at supersonic points, these terms 
represent a positive artificial viscosity which, when the flow Is 
aligned with the X direction, reduces to 



c^) 4X 


^XX 


as In the scheme of Murman and Cole [20]. 

Equations (38) to (42) are used at each Iteration of a run. To 
describe the iteration process, which Is done by row relaxation, it Is 
helpful to Introduce an artificial time t, which Increases by the 
quantity At at each Iteration. The right-hand sides of equations (38), 
(39) and (41) become 


®l+l,J,k" 


2(l-l/u) - (2/dj) 


(4X)' 



f 



0 

1+1, j+l,k 


pO . pN 

®l-l,j+l,k -<5i+l,j-l,k+ Gl-l,j-l,k 
4 ax ay ” 



2G 


N 

i.j»k 


rO 

®l,J,h 


2gJ. 


l.J.h 


+ G 


0 

i-2.j,k 


(AX) 2 



vAiere the superscripts N and 0 denote new or old values of the 
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potential and where w Is the overrelaxatlon factor. These expressions 
represent approximations to 


"XX 


At 

AX 


^°Xt 


AX u 


OG^). 


'XY 


I At ,, 


"XX 


+ 2 


~ Gv. 

AX 


rfence the equation solved Is 

(43) I (c2p^.- BiBOGXiX. + I 0^)0^^+ a^Gx^J +60^-0 , 

i,j ^ 


where the coefficients and 5 are determined by the new 

approximations (38), (39) and (41). 

By considering an orthonormal system of coordinates (p,r,s), where 
the flow direction is still the s direction, we arrive at the equation 


(4) (c2-q2)G5^ + q^r^Opp + 5iGst+ 82 G„+ 6,0^^ + 6G,. = 0 


If 


" ’ ' 2 ® +-T ■'“2 

^ c^-q^ q^ q^ 


is a new time coordinate, then equation (44) is transformed into 


(45) (c2-q2)G 


ss 


+ q2c 


rr 


9 ®pp - (• 


01 


2 2 
c^-q"^ 




-4>eiT 

q“ 


+ 6G, 



In order to ensure the convergence of the iteration scheme to a 
solution of the steady-state equation, we want equation (45) to be a 
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daniped three-dimensional wave equation. Given the form of equation 
(45) this Is equivalent to the condition 



At subsonic points the damping condition is always satisfied. At 
supersonic points it may no longer be true. The choice of 
approximations (41) and (42) help to ensure a large coefficient 0^ , 
but near the sonic line this may not be enough. One way to ensure that 
condition (46) Is satisfied, Is to Increase 0 ^ by adding a term of the 
form 


0 — 0 . 
AX St 


® ^ ^ M ®2^Yt‘*' ®3^Zt^ 


to equation (43), where 0 Is a positive parameter large enough to 
•Increase 0 j by the right quantity. 

The mixed derivatives used above, Gj^^. , Gy^ , are obtained by 
using approximations of the type 

pN pO «N 4 . pO 

p . ®i,j,k “ ®l,j,k Gl-l,j,k + Gi-l,j,k 

Gxt 3 TTx 


Me have described the difference scheme for Interior points. For 
the boundary conditions at the hub and cowling, we transform equation 
(35) Into a difference equation In order to compute the value of the 
reduced potential at these boundaries. For the blade surface and the 
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vortex slieet we introduce ghost points behind tfie boundary. These 
points are used to compute the value of the potential at points on the 
blade as if they were interior points. The values of the potential at 
the ghost points behind t)« blade are determined by using equation 
(36), i.e. the slip condition. For the points behind the vortex sheet 
the potential is determined by the Kutta-Joukowskl condition. To 
compute the value of G at the remaining boundary points we use an 
approximation to the outlet velocity. 

The program permits not only compressor blades but also propellers 
to be modeled. Since propellers have no cowling, the inlet speed is 
equal to the outlet speed. Thus, for a propeller run, all the boundary 
conditions except those at the hub, on the blade and on the vortex 
sheet are replaced by 

G = 0 . 

With program FL022 as the starting-point, we have used the result 
stated in this section to write a new program in FORTRAN named CSCDF22. 
This program Is listed in Section VIII. 
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V. RESULTS 

In the preceding sections we described our computational method. 
In order to show how this method wocRs, we present in this section 
seven exanples of runs made with the program CSCDF22. 

The first case was run to test the validity of our code. A 
compressor blade had been designed with code R [3]. Fl}»ure 2 shows the 
result of this run. We used the coordinates of this two-dimensional 
blade profile to create a three-dimensional blade. The profiles were 
identical at each span station. Since we chose to run the program for 
a cascade configuration and in the compressor case, we were solving a 
two-dimensional cascade problem. The Inlet speed and the Inlet angle 
were the same as In the design run, but the blade was markedly cambered 
and we chose a gap-to-chord ratio of 1.5. Figures 3 and 4 show a 
representation of the cascade and the grid at the hub which Is the same 
at each span station. The pressure distribution Is given In Figure 5. 

The other six examples are all run with the same blade. The 

profile Is defined at three different span stations. The bosic profile 
is the NACA-0012 profile. At each span station the chord and thickness 
are the same but, for the realistic axial flow configuration, the 
coordinates of the singular points are different. The blade has a 
sweep angle of 14 degrees and the dihedral angle decreases from 0. 
degrees at the hub to -10. degrees at the tip. The angle of twist is 

equal to 0. degrees at the hub, 2. degrees at the cowling, and 1. 

degree in between. For the axial flow configuration the distance 
between the hub and the axis Is 2. This Is also the distance between 
the hub and the cowling (or the tip). 

The first example with this blade Is a run for a compressor 
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cascade configuration which appears In Figure 6. Figures 7 and 8 show a 
representation of the grid and the pressure distribution. The 
different flow parameters are shown In Figure 8. Ml is the inlet Mach 
niunber. M2 Is the outline Mach number. DEV Is the difference between 
the outlet angle and the inlet angle given by ALP. 

The four next runs are also analyses of compressor flow, but now 

In an axial flow configuration and with eight blades around the hub. 

For each run the parameter OM determining the speed of rotation has 

different values. These values are 0.0, 0,5, l.O and 2.0. The other 

flow parameters are identical. Figures 9 and 10 represent the grid at 

the hub and cowling. It will be noticed how different they are, 

although the profile Is almost the same at the hub and cowling) but the 

gap-to-chord ratio is twice as large at the cowling. Figure 11 shows 

how the blade looks in a plane orthogonal to the axis. The different 

pressure distributions obtained for each run are given In Figure 12, 

13, 14 and 13. These figures clearly illustrate how the flow evolves 

as the speed of rotation Increases: the shock on the upper surface 

/ 

weakens and then disappears while, on the other hand, a shock appears 
on the lower surface and intensifies. 

The last example is similar to the fifth, except that it is a 
propeller analysis. Figure 16, 17, and 18 show the geometry of the 
case, and the pressure distribution Is given in Figure 19. by 
comparing Figures 14 and 19, we observe that, for example 5, the shock 
appearing at the cowling which Is caused by the speed of rotation. Is 
amplified significantly by the cowling. 

These numerical experiments show that our computational method, in 
spite of Its restrictive hypotheses, enables us to analyze the 
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thcee-di:nenB tonal transonic flow around compressor blades. This Is a 
first step towards the study of compressors or propellers in three 
dimensions. Considering what happened In the study of transonic flow 
past swept wings, we venture to suggest that the next steps could be to 
improve this method In order to Include the computation of the wave 
drag and to design blades with a prescribed pressure distribution 
[4. HI. 
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VI. HOW TO USE THE CODE 

This section Is Intended to serve as a guide for users of our 
computer code. We explain the different code options and give a 
listing of the Input parameters. 

Our computer code Includes two options which permit the solution 
of either the compressor problem or the propeller problem, and It 
treats them either In a realistic axial flow configuration or In a 
three-dimensional cascade representation. In the first option there Is 
a choice of two boundary conditions. In the propeller case the reduced 
velocity potential Is set to zero at Infinity; In the compressor case 
the boundary condition Is given by a Neumann condition at the cowling. 
The first option determines the first mapping to be performed. For the 
axial flow configuration the whole space Is mapped Into a part of the 
space consisting of as many periodic strips as there are blades around 
the hub. For the cascade configuration this first mapping Is not 
required, since the computation starts directly from a periodic strip. 
In both cases the periodic strip Is mapped conformally Into a 
quadrilateral. For the compressor case the blade shape and the vortex 
sheet define the bottom of the quadrilateral; the hub and the cowling 
determine the sides; and the upper edge Is defined by the flow at 
Infinity. For the propeller case the side previously determined by the 
cowling Is defined by the flow at Infinity. 

The shape of the blade Is defined by giving the coordinates at 
different span sections from the hub to the cowling. As many as six 
sections can be defined by giving the Cartesian coordinates at each 
section with z constant. If two sections are similar, only the 
coordinates of the first section are read In. . 
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W* shall now describe the input dattt cards and explain the 
paraiaetecs which occur In the input. All the data cards are read In 
fomat (8F10.6). A title card is insert«iid before each data card or 
each aeries of data cards. 

Title Card 1: 

NX, NY, NZ, FPLOT, XSCAL, PSCM., GRIDX, GRI0Y. 

NX, NY, NZ are the number of mesh cells In each direction of the 
transformed spec#. Thus NZ is the number of mesh cells 
In the radial direction, NX the number of mesh cells 
along the blade and the vortex sheet, and NY the number 
of mesh cells in the third direction of the 

computational grid. If NX ■ 0^ the program stops. The 
dimensions in the code are NX+l , NY+2, NZ+3. 

FPLOT is the parameter controlling the generation of plots. 

IF FFLOT ■ 0, a printer plot but no Calcomp plot is 
obtained at each span station. If FPLOT • I, both a 
printer plot and a Calcomp plot are generated; and if 
FPLOT - 2, only a Calcomp plot is generated. 

XSCAL, PSCAL determine the scales of the Calcomp plot. The pressure 
scale is set to PSCAL per inch in each section plot. 
PSCAL ■ 0 is equivalent to PSCAL "0.5. If PSCAL is 

positive, each section is scaled so that the span is 5. 
If PSCAL is negative, each section is scaled 
proportionately to the local chord and the maximum chord 
is XSCAL/2. 

GRIDX, GRIEV control the grid generation. The grid depends largely 
on the geometry of the compressor, so these two 
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Tltle Card 2 
MIT, COV, PI 
MIT 
COV 


P1,P2 

P3 

SIA 

FHALF 


paramecers nuat be chosen carefully so as to obtain an 
acceptable grid. GRIDX has to be greater than -I and 
GRICflf Is positive. In most of the cases we can use 
GRIDX - 0 and GRIDY - 1. Decreasing GRIDX gives more 

points near the nose and less at the trailing edge. 
Increasing GRICK gives fewer points near the nose and 
more at the trailing edge. Increasing GRIDY Increases 
the number of points near the blade and the vortex 
sheet . 

. P2, P3, BETA, FHALF, FCONT. 

is the maximum number of iterations computed, 
is the desired accuracy. If the maximum correction is 
less than COV, the program goes to the next grid or 
terminates. The same will happen If the .maximum 
correction is greater than 10. 

are the relaxation parameters for the subsonic and the 
supersonic points respectively PI Is between I. and 2., 
P2 Is less than or equal to 1. 

Is used for the boundary values. P3 Is usually set to 
1, but a smaller value Improves the convergence for very 
cambered blades. 

Is the damping parameter controlling the amount of added 
♦ gj.. It is normally set between 0. and 0.5. 

Indicates whether the mesh Is to be refined. For FHALF 
■ 0 the program stops after convergence or after the 
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FCONT 


Title 

FMACH 

FMACH 

OM 


AL 

PA 

ZONE 


prescribed number of Iterations have been completed. 
For FHALF - 1 the mesh will be refined and another input 
card must be read for the parameters of the refined 
mesh. The mesh can be refined only twice and the 
maximum grid Is 24x14. 

Indicates how the computation starts. After each run 
the velocity potential is stored on tape 8. If FCONT ■ 0 
the program begins by initializing the velocity 

potential. If FCDNT « I , we have a continuation run and 
the data stored on tape 8 from a previous run Is used as 
input and is read in on tape 7. 

Card 3: 

OM, AL, PA, ZONE, DC. 

is the Inlet Mach number. 

is the speed of rotation. This means that the 

compressor is turning at the speed of OM revolutions per 
60 units of time. The unit of time depends on the flow 
conditions. The unit of length is the length of the 
blade divided by CHCXID and the unit of speed is the 
Inlet speed. 

Is the inlet angle. 

determines whether an axial flow configuration or a 
cascade is to be studied. PA ■ 0 for a cascade. PA ■ 1 
for an axial flow configuration. 

Is the distance between the axis and the hub. In the 
cascade case It Is set equal to 1. 
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DC Is the gap-tp-chord ratio for the cascade case. For an 

axial flow configuration it is the number of blades 
around the hub. 

Title Card 4: 

ZSYM, NC, SWEEP!, SWEEP2, SWEEP, DIHEDl, DIHED2, DIHED 

ZSYM indicates whether we have a compressor problem or a 

propeller problem. ZSYM ■ 0 for the compressor case and 
ZSYM ■ 1 for the propeller case. Thus, if ZSYM ■ I, 
there is no cowling, the outlet speed is the same as the 
inlet speed and the radial direction is stretched to 
Infinity. 

NC is the number of span stations at which the blade 

section is defined. The blade section is interpolated 
linearly between two span stations. If NC < 3 the 

geometry is given by the last case studied and the 
program will run for this ZSYM and for the parameters 
given by data cardc. I through 3 without the need to read 
any more cards. 

SWEEPl , SWEEP 2, SWEEP are the angles of sweep of the singular line at the 
hub, at the cowling (or the tip) and at the far field 
respectively. For a compressor blade SWEEP is not used. 

DIHEDl, OIHEKt , DIHED are the dihedral angles of the singular line, at 
the hub, at the cowling (or the tip) and at the far 
field respectively. In the compressor case DIHED is not 


used 
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Title Card 5: 
Z, XLE, YLE, 

Z 

XLE,YLE 

CHORD 

THICK 

ALPHA 

FSEC 


Title Card 6: 
YSYM» NU, NL 
YSYM 


NU,NL 

Title Card 7 


CHORD, THICK, ALPHA, FSEC 

Is the location of the span section, 
are the coordinates of the leading edge. 

Is the chord of this section, used to scale the profile, 
nultlplles each y coordinate; thus the thickness of the 
sect ion is nul t ip lied by THICK. 

is the angle of twist through which the section is 
rotated. Changing the value of ALPHA by the same amount 
at each span station Introduces a stagger angle. 
Indicates whether section coordinates are to be read In 
from data cards. FSEC ■ 0 means no further cards are to 
be read for this section. FSEC ■ I means the 
coordinates of the section will be determined by the 
next Input cards. 


indicates whether the profile Is symmetric or not. YSYM 
■ 1 for a symmetric profile and YSYM ■ 0 otherwise. For 
YSYM ■ 1 only the coordinates of the upper surface are 
read. 

are the number of points on the upper and lower surfaces 
respectively. 


TRAIL, SLOPT, XSINB, YSINS 
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TRAIL 

SLOPT 


XSING 


Title 
X, Y 
X.Y 


Is the Included trailing edge angle In degrees, or the 
angle between ^he upper and lower surfaces if the 
profile is open. 

is the slope of the mean camber line at the trailing 
edge. SLOPT is used not only for the profile but also 
to determine the vortex sheet behind the profile. If 
SLOPT is too big and the gap-to-chord ratio is small the 
program may abort. 

YSING are the coordinates of the singular point inside the 
profile about which the profile is unwrapped. This 
point nust be chosen in such a way tliat the mapped 
coordinates are as smooth as possible without a 
pronounced bump near the nose. If such a bump occurs, 
we first move the bump nearer the nose by positioning 
the singular point closer to the upper or lower surface. 
Then the bump is smoothed by placing the singular point 
at the right distance from the nose. Moreover, the 
singular point must be redefined not only for any new 
profile but also for any new configuration, since the 
geometry changes with the parameter DC. 

Card 8: 


are the coordinates of the upper surface. There are NU 
cards and the coordinates are read from the leading edge 
to the trailing edge. 
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Tltle Card 9: 

X, Y 

X,Y are the coordinates of the lower surface. There are NL 

cards and the coordinates are read again from the 
leading edge to the trailing edge. The first data card 
repeats the first data card for the upper surface which 
defines the coordinates of the leading edge. 


Ihese cards are followed by (NC-1) series of cards beginning with 
Title Card 5. 

The program provides both graphical and printed output. The 
Calcorap plots were shown In the previous section. 

In the printed output we can read the listing of all the 
coordinates given In the Input. Then the mapped surface coordinates 
are printed for the hub and the tip. This allows us to ascertain 
whether the coordinates of the singular points are well chosen. After 
the iteration history, with the maximum correction and the mean 
correction to the velocity potential and with the maximum residual and 
the mean residual for the difference equations at each iteration, the 
program gives us a printed plot of the coefficient of pressure on the 
surface. This is followed by a Mach number chart. These last two are 
given for each span station. Finally the characteristics of the blade 
are printed . 

If the mesh has to be refined, a new series of output is printed 
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for the new raesh. If not, the program terminates or restarts by 
reading a new data deck. 
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original PAGE m 

OF POOR QUALITY 


VIII. LISTING OF THE CCOt 


An CSCDF2 2( INPUT* 512 » OUTPUT- 51 2> TAP 6 6 ■INI'UT# TAP E6-UUTPOT# TAPE 
17*5 J,2» TAfE5»5l2) 

C NAIN PjuTINE controls TH6 COM PUT A T I ON A L PROCEDURE 

C G IS He »• EDUCED VELOCITV POTENTIAL 

C'ln'^lN b(l29»2t,17)»SOa29»17)#EO(17>» IV(129»17HITE1(17) »ITE2(I7) 
l,AO{129)>Am29),A2(i2 7)^A3(129)»BO(26),ail26)»B2(26)»fl3(2Dr>Z(i7) 
? , C 1 a 7 ) # C2 ( 1 7 ) » Ci(17 ) » C <1 7 )> C 5 ( 17 ) » KC ( i 7 ) . X 2 (17 ) , X2 Z ( 17 ) , YC ( 1 7 ) » f 
'iZ ( 17 ) • yz 2 (17)»KSYH»NX»Ny»N2#KT6l»KTE2^IiyMiSCAL»SCALZ# XK » OnE G A » A L P 
AHA.CA, SA»MACH»T>PA, ABLAOE 

common 7CAL/ Pl»P2;P3>3ETA»pR^K'1»DG#Gn»NS.Ul»Vl#Wl/IR»jK>Kk» I3)*JG» 
IKG 

QIMENSIJN XS(20J#6)» YS(200#6»# ZS(6)» XLt(o)# YLE(6)> SLQPT(6>* T 
1KAIL(6)» NPUH E1(6)» E2{6)» Ei(6)» EA(o)» E5(oH XP( 200)» YP{20C 
Z)f JJ(200)» D2(2CO)» C3( 200)» X(129)# Y(129)f SV(129)» SNU29), C® 
3(12-?)# !.HDRU(i7)> SCL(17)» SCD(17)# SCn(17)» PITH)# CQVO(3), P1L( 
A3). P20(3)» PiC(2)» 36TAC{3)» FHALF(3)» RES(501)> C0UNT(60i)» DE SC 
5 n 0 ) » n T L E ( i 0 ) 

NO»2uu 

NE*129 

INu»l 

Ul*C. 

Vl-0. 

». 1 * 0 . ' 

HE AO = 3 
IhRIT» 6 
KPLOT*J 
IPLOT»l 
J I T * 0 
(JEwlNO 7 
RAD*£'7. 29578 
10 WRITE (HkIT.430) 

WRITE (URIT.210) 

READ (UEADH20) TITLE 
WRITE (IWRITH60) TITLE 
READ (iKEAD»A20) DESC 

READ (iKtADHio) fnx,fny,fnz.fplot.xscal.ps:al.f»ty 
IF (FNX.lt. 1. ) GG TO 200 

WRITE (IwRIT.620) FN X . FN Y . FNZ » F PL OT # X SC AL » PS C AL > F» T Y 

NX-FNX 

NY«FNY 

NZ = FiNZ 

KPL0T*ABS(FPL0T) 

READ ( IREA0.A20) DESC 
WRITE ( IWR IT. 530) 

NM*0 

20 NM*NM + 1 


ORIQINAL PAGE 
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READ < IRE AO# ^10) F I T< N« ) #C 0 VO ( Nfl) # R 10 ( NM ) » P 2 0 (NM ) # P 30 ) # d E T A C-( N« 

1 ) »FHALF<Nrt>,FCONT 

WRITE (UP IT #‘♦50) FIT(NM)#CUV0(NM)#PlC(NM)#P20(NM),P30(NK)#dETA0(N 
IMUFHALFiNM) 

IF (FHALF(.SM).NE.O.. ANO.NM.lt. 3) GO TO 20 
FHALF(3)-0. 

READ (1REA0#A20) OESC 

READ (IR6A0#A1C) F MAC H# DM# AL » P A # ZON 6 # DC 
WRITE (UPIT»5AC) FM AC M, OM # AL # P A » DC 
DC* (X. -PA) ♦DC ♦PAD/ 360. ♦PA/ DC 
XKwl./OC 

ZDNE'i.^PAW ( ZONE-i. ) 

QMEGA«PA*aM/RAD+6. /ZONE 
ALPHA*AL/R AO 

CALL G6QM (ND#NC#NP#ZS#)(S# TS#XL6# VLE# SLOP T# TR A IL # X P# YP# S W E tP 1 # SwEE 
1P2# SWtEP#0 1 HED I# DIHt02# 01 HEO#KSYM#XTEO# CHQRDO#ZT IP# I SY‘1C# XK# PA# INC 

2 ) 

ISYM*U‘YMO 

IF ( alpha, NE .0. ) ISYM-0 
CA«COS (ALPHA ) 

SA«SIN(ALPHA ) 

IF (FCGNT.LT.1) go TO 40 

READ (7) NX#NY#NZ»NM#K1#K2#JIT#U1»V1#WI 

MX-NX+i 

MY-NY + 2 

MZ-NZ+i 

DO 30 K-1#MZ 

READ (7) ( (G(I#J#K)#I«i#MX)#J*l#MY) 

30 CONTINUE 

READ (7) ( E0(i<)#K»Kl#K2) 

40 CONTINUE 

CALL COORO (nx#ny»nz#xteo#ztip#xmax#zmax#ksym#pa#sy#scal#scalz#ax# 
lAY#AZ#Ty#F#AO#Al#A2#A3#0O#Bl#a2#B3#Z#Ci#C2#C3#C4#C5) 

CALL SINGL (NC#NZ#KTE1#KTE2#CH0RD0# PA# SWEEP!# sweep 2# sweep #D lht01#0 
1IHED2#0IHE0# ZS »XLE#YLE#XC#XZ# XZZ# YC# YZ» YZZ#Z#6l#E2#E3»E4#E5# INO»ZO 
2NE ) 

CALL SURF (ND#NE#NC#NX#NZ# ISYM#KT£1#KTE2# SCAL#AO#Z#ZS#XC# YC#SLO‘’T# 
ITRAIL#XS#YS#NP# ITE 1# I TE2 # I V# SO# XP # YP# Di# 02# 03 # X, Y# IND#XK#PA# XZ» YZ# 
2Ail#Cl#KSYM) 

IF (INO.EO.O) go to 190 
IF (FCONT.GE.l .) GO TO 50 
NM-1 

CALL ESTIM 
50 CONTINUE 
NIT -JIT 
T-2./SCAL 

60 WRITE (IWRIT#430) 

FCONT-0. 

MiIT*FIT(NM)+JIT 

COV*COVO(NM) 

BFTA«BETA0(NN) 

MX-NX+1 


«y 

MZ«N 2*3 

KY*NYfi 

Kl*3 


*»»IU (lW^IT»/i20) 
^^•'] 7C I-2#NX 

(UKir^^eo) 

7C CaNflNUc 


uoric (X*FIT»43C) 
">% I T t ( I 1 T # 2 30 ) 
C3 6u 


w^ITb ( Uf I F 
PO CDNTlNUb 

h-iiT-: (UPir 


i*>. i T h 
Mi TF 
i.Kr rf. 

DO 90 

re 


( 1 w M T 

( UF IT 
r< y 
( UF £T 


440) 

240 
440) 
4 30 
250) 

4 40) 


ORIGINAL PImd ^ 
OF POOR QUAUTV 


nv< I »<) #K»K1 »K2) 


A0( I )»S0(lrKl)»SuU»KTE2) 


XHAXMX 


80 ( J) 


•/C 


cinh sut 

.♦iPITfc (UK II 
WPI Tt 
W^^ITE 
WiPITt 
OQ 100 
M'i Tfc 
wP.ITb 


( UP IT 
( UP i T 
( I .*kIT 
K, » K U K 
( UP IT 
(U-^IT 


100 


C 0 NTINO 6 
WPITE nWFIT 


( UKIT 
(UR IT 
( UP IT 
(UKIT 
( UR IT 
SECOND 
( IwF IT 
( UP IT 
( URIT 
(UKIT 


260 
440) 
430 ) 
27C ) 
2 

440 ) 
440 ) 


SYf AY 


Z (K)»XC(K).YC(K)#XZ( K) »Y2( K)»XZ/(K),YZZ(x) 


280 ) 
440 ) 

4 30 ) 

2 90) 

3 00 ) 
470 

(TIME) 

5 10 ) 
310) 
440 
320 ) 


ZMAX» AZ 


NX»NY»NZ 


TIKE 

FMACH»OM» 4L 


110 


WRITE 
WR I I E 
MITE 
wRI TE 
WRITE 
CALL 
K R I T E 
WRITE 
WRITE 
WRI TE 
LX*NX/2+l 
CL*0. 

DO 110 K*Ki#K2 
Il-ITEKK) 

X( ID »U+.5tSCAL*( A0( I1)YA0( Ii)-SO( I1,K )*S0( rUK ) ) 
YUSCAL*A0(I1)«S0(IUK) 

X (I1)»XC (K ) + ALOG(K (Il)*W2 + YU«2)/XK/2. 

X(LX)*1 .♦.3»$CAL^( AO(LXH‘AO(LX )-S0( LX>K)^SO( LX#K ) ) 
YUSC4L*A0 (LX)«SC (LXrK ) 

X (LX ) -XC (K) ♦ALOGC X ( LX )**2>YU*2 )/XK/2 . 

CHQPt(K)*X( Il)-X(LX) 

CONTINUE 


ORIQIMal msa 
OF FOOf? QUAUry 


KZ0U^»KTb2-l 

S>0. 

DO 120 K«Krei»KZ0UM 
0Z0«.;?M/(K + 1)-Z(K) ) 

S-S + CZC^ (C40R0<Kfl)4*CHOi?0(K) )/C5(K) 

120 CCNTlNUb 
ABLADE-S 
130 NIT-NIT+1 
P1»P10(NM) 

P2»P20(NM) 

P3*P30(NH) 

CALL BOUND 

12-Ui+CA 

V2-VUSA 

WRITE ( IWPIT»^90) N1‘T»0-3»IG# JG#KG»FR» IR> JR»<R#GM#«,'i>FlfU2»V2^NS 
COUNTiNIT) »NIT-1-JIT 
RES (Nin -FR 

IF (NlT.Lr.Mir.AND.ABS<0G).GT.C0V.AND.AB5(DG).LT.lC.) GO TO UO 
RATE«0. 

IF (NIT.GT.JIT + n RATE»{ABS(RES(Nm/RES(JIT4.1)))ft( 1 . / ( C JUN T ( N I T ) 
1-C3UNT (J lT + 1) ) ) 

WRITE (I,<Rir# 330) 

WRITE (IWRIT» 500) RE S ( J I T^- 1 ) # RE S C NI T ) # C OUNTI N I T ) , R ATE 
CALL SECOND (TIME) 

WRITE (IWKITolO) TIME 
LX-NX/ 2+1 
DC 150 K»K1.»MZ 

IF (K.L T.KTE1.0R.K.GT.KTE2 ) GO TO 150 

ll-ITEKK) 

r2,*lTE2<K) 

CALL VELD (KiK, SV»SM,CP#X, Y) 

CH0RD(K)-X(I1)-X(LX) 

CALL FORCE ( 1 1 > I 2 # X » Y # C P # A L # CHQ RO ( K ) » XC ( K ) » S C L ( K ) , SC D U ) * SCM (K ) # A2 
1, A3,C5(iO#PA) 

IF^ (KPLOT.GT.l.AND.K.GT.KTEl) GO TO 1<»0 
WRITE (1WRIT#A30) 

WRITE (IWRIT#3A0) 

WRITE (IWRIT»AAO) FMACH»OM»AL 
WRITE ( IWP IT#350) 

lAO WRITE (IWRIT»A40) Z (K )» SCL (K )» SCO ( K )# SCM( K ). CHOP D( K ) 

IF (KPLOT.LE.l.) CALL CPLOT ( 1 1# 12, SM» AO, S0( 1 #K ) , FMACH) 

CALL SPEED (K) 

150 CONTINUE 

CALL TOT FDR ( KTE1,KTE2,CH0RD, SCL,SCD, SCM, Z,XC,C5,CL,CD,CMP,CMP,CMY 
1,PA,A3LADE) 

VLD-0. 

IF (ABS(CD) .GT .l.E-6) VL0-C4./C0 
WRITE (IWPIT,A30) 

WRITE (lwRIT,360) 

IF (KSYM.EQ.U U2«CA 
IF (KSYM.EQ.l) V2«SA 
V«ATAN2(V2,U2) 
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OF PC- 


.u ..i 

:auw 


V«V*CAu-AL 

u*sopr( w/( 

kPITc (Iv.Pn>AA0) fmach»oh»al>u# V 
(Ik?IT»370 ) 

kPITE (1,^5 IT# A AO) CL»C0,VLD 
kRlTt (1 nPIT»380) 

kRITfe (£taPIT#AAC) C«P#CM>^#CMY# ABLAD6 
LPL0T.pPLOT^( I .-FH4LF(N"1) ) 

IF (LPLOT.lt. 1) GQ TO L60 * ■ 

DC«360./XH/PAD*ZCNF/CH0k00 

CALL TH-^cEO { IPLOT»SV»S-'i#CP#X#y»TITLE»DC# AL# ZC!NE#U#V#CHDPOO# X,tCAL» 
IPSC AL) 

160 CONTINUE 

IF (FhAlF( '(M l.fcC.C.) GO TO 170 
NX« NX4 -n;^ 

ny-ny+ky 

NZ»NZfNi 

CALL CCOPD ( NX»NY#NZ# XT60» ZTIP»XMAX#ZMAX»KSYM»PA#S Y> SCAL# SCALZ# AX# 
lAY»AZ»TY#FfA0»Al»A2»A3#B0#Bl»32»d3#Z#Cl#C2»C3»CA»C5) 

CALL iINGL (NC#NZ#KTEl#KTE2#CHQi?00#PA#Sw£EPi#SW££P2»SWfc£p»DIHE01#r/ 
lIHeC2#CIH60# ZS»XLE# YLE# XC# XZ# XZZ» YC#YZ#YZ2»Z#61#E2#E3#EA#E5» INO# 20 
2NE) 

CALL SUPF (NDrN6#NC#NX.NZ» ISYM»KTE1#KTE2# SC AL * AO i Z # Z S # XC # YC# SL OP T# 
ITRA IL# XS»YS# N?» ITEl# ITE2# IV#S0»XP#YP#01#02#0 3#X# Y# iNOrXK^PA^XZ^YZ# 
2A1#C1#KSY^) 

IF (IND.EC.C) GC TG 190 
CALL REFIN 
NM.NM+1 
N1T»C 
GO rc 60 
170 CONTINUE 

WRITE (6) SXfNYfNZiN'I^KUKZjNlTfUl# VI# W1 
WRITE (6#390) 

DO leo K«i#Mz 

WRITE (8) { (G( I# J#K) » I«1»MX)# J -ItNY ) 

180 CONTINUE 

WRITE (8) (EQ(K),K*K1#K2) 

END FILE 3 
REWIND 0 
GO TO 10 

190 WRITE (1WRIT»930) 

WRITf (IWRIT,400) 

STOP 0102 
200 CONTINUE 

CALL PLOT (0.#0.#999) 

STOP 0101 

210 FORMAT (20H0NYU COMPRESSOR CODE) 

220 FORMAT t AONOINOICATION OF LOCATION OF BLADE AND VORTEX SHEET#274 1 
IN COORDINATE PLANE Y • 0. / 27H0 ( UV ( I#K ) #K«K1 # K2 ) # I *2 # NX ) ) 



ORIGINAL lU 
OF POOR QUALIW 
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230 fOkMAT <‘»9H0CHDR0wIS£ CELL OISTRIdLTIDN IN TRANSfQkMEO PLANE, A 
INO HAPPtC SURFACE COORDINATES AT HU8 AND TI^/J5H0 > ,1^ 

2H HUB PHQFILE>15H TIP PROFILE ) 

2A0 FORMAT (16HG TE LOCATION ,i5H POWER LAW ) 

260 FORMAT (WMhONDFMAL CELL DISTRIBUTION IN TRANSFORMEC PLANE/16H0 
IT) 

260 FORMAT (liJHO SCALE FACT0R»15H POWER LAW ) 

270 FORMAT ( W^HOS^ANWl SE CELL DISTRIBUTION AND SINGULAR LINE/16HC 
I I #I6H X SING ,16H Y SING , 15H XZ 

2»16H YZ ,16H XZZ ,15H Y7Z ) 

200 FORMAT (I6H0 TIP LOCATION, 15H POWER LAW ) 

290 FORMAT ( i9H0ITERA TI VE SOLUTION) 

300 FORMAT (I6H0 NX ,I5H NY ,15H NZ ) ‘ 

310 FORMAT (IpHO MACH NO ,I5H OMEGA ,15H ANG OF ATTACK) 

32C format ( lOHOITERAT ion, 16H CORRECTION , AH I , AH J ,AH K ,15h 
I RESIDUAL »AH I ,AH J ,AH K ,10H MEAN C0R.,10H MEAN RES., 10 
2H REL FCT XSPEED , ICH YSPEED ,10H SONIC PTS) 

330 FORMAT (15H0 MAX RESIDAL 1,I5H MAX RESIQAL 2»15H wQRK ,i ) 

16H PEDUCTN/CYCLE ) * 

3A0 FORMAT (2AHCSECTI0N CHAR AC TE R I S T IC S / 1 5H0 MACH NO ,15H OM ’ 

lEGA ,15H ANG OF ATTACK) 

350 FORMAT (/13H SPAN STAT ION, 12X2HCL, 1 3X 2HC0, 13X2HCM, 10X5HCHQRD ) 

360 FORMAT (22H03LA0E CH AR AC TER I ST I C S/ 1 5H0 MACH NO 1 ,15H QMEG > 

lA ,I5H ANG OF ATTACKflSH MACH NO 2 »15H DEV. ANGLE ) 

370 FORMAT (15HC CL ,I5H CD ,15H L/D ) 

36C FORMAT (/ZXFHCM PITCH, 6X7HCM R0LL#9X6HCM Y Aw , 9X 6HABL ADE ) * 

390 FORMAT (1X»1AHWRITE ON TAPES) ( 

400 FORMAT (24H0BA0 DATA, SPLINE FAILURE) 

410 FORMAT ieFiO.6) I 

420 FORMAT (1OA0) 

430 FORMAT (IHl) I 

440 FORMAT ( FI2 . 5 , 7F1 5 . 5) j 

45C FORMAT (IX,8E15.5) I 

460 FORM/T (1HO»10A3) 

470 FORMAT (10,7115) S 

400 FORMAT (IX, 3214) | 

490 format ( 110,E15. 5,314, E15. 5,314, 2E10. 3, 3F10. 5,110) 

50C FORMAT (2E15.4,2F15.4 ) | 

51C FORMAT ( 15H0C JMPUT ING TIME , FIO . 3, lOH SECONDS) 

520 FORMAT ( 75 X3HFN X, 1 IX 3HFNY, 11 X 3HFNZ, lOX 5HF PLOT ,9X 5HXSC AL, 9X5HPSC AL , S 

19X5HGP1DX>9X5HGRIDY/1X,8E14.5) # 

53C FORMAT ( /4X7HF IT < NM ) , 8 X8HC0Y0 ( NM ) # 8 X7HP10 (NM ) ,8 X 7HP20 (NM ) , 0 X 7HP 30 ( - rf 

lNM),6X9H6ETA0(NM),6X9hFHALF(NM)) j 

540 FORMAT (/5X5HFMACH,llX2HaM,14X2HAL# 12X3H PA, 12X4H DC /iX,5E15.5) I 

END I 

ii 

? 


I 


I 

i 

I 
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OF POOR i 


SUBI<CUTIn& GBQfl (NO/NC . ZS, X S# YS » XL £ # U £ # i L GP f » T« A I L » X P, YP,i».rEP 
ll#S«ecP2#i^»:6P>DIHE01»DIHPD2^0IHED#KSYM»XTE0»CH0PU0, ZT I F , IS YMO » X K, 
2PA,INO) 

GEJAfCr-^IC OEFIMriON CF SLAUE 

DIMfcNSIJN KS(^^D^1)# YS(N0.1)> ZS(1)> XLE<1)» YLE(1)» SluPHUf ?T( 
16)» 1RAU(1)» XP(1), YP<1), NP(D# K(6)» Y(6)» Z(6)» XC(b), YC(t)» 
2 C(6) 

DIMENSION OtSC(lO)/ A(o), d(t), D(6), E(6) 

IREAD*;? 

IWRIT»6 

8AD-57 .2<;i? Zfl 

READ (IR£Aa#23C) DESC 

READ ( IkEAO# 22C ) 2SYM,FNC# SwEEPi# SWEEP 2» 5WEE P »DIHEDl»0lHED2,0I hF.C 
taRITF ( I WRIT# 2 70) ZS YM » FNC ^ S wEE P 1# S WEE P2» S WE E P# 0 I HF 0 1 # 0 1 HE 02 » D H E D 
RS Y^'^SYM 

IF {FNC.LT.3.) RETURN 
NC-FNC 

SwEERl»SwtEPl/RAD 
SWEEP2>SWEEP2/RAP 
SWEEP»iwEEP/RAD 
DlHEDi=DIHEDi/ RAD 
0IHEC2*D1HEU2/ PAD 
DIHEC-DIHEO/RAO 
S1»TAN (SWFEPI) 

S2*TAN( S»rEEP2) 

T1»TAN(DIHE01) 

T2»TAN (DIHED2) 

ISYrC»l 

XTEC-0. 

CHORD C=0. 

ZCNE-1. 

K-1 

1C READ (IREAD#230) DESC 

READ (IREAD>220) Z T ( K ) ^ X L# YL > CHORD# TH I C K» AL » r SEC 
WRITE UWRIT#280) Z T ( K )» XL # YL# CHORD # TH ICK , AL # FSEC 
alpha»al/rad 

ZriK )«ZT (K)/ZQN£ 

ZS(K)»aLOG(ZT( K) + (1.-PA) IwPA+d.-PAIWZnO 
IF (FSEC.£a*0« ) 60 TO 80 
IF (K.GT.l) GO TO 20 
Z0NE«(1.-PA)+PAWZT(1) 

ZS( 1 )*( I .-PA) WZII 1) 

ZT(1)«1. 

20 CONTINUE 

READ ( IRfcA0#230 ' DESC 

READ (IREA0>22C) YSYM#FNU#FNL 

WRITE UWPIT#290 ) Y S YM# FNU# F NL # ZONE 

NU-FNU 

NL-FNL 

IF ( YSYM.EQ.l* ) NL»NU 
N«NU+NL-1 

READ (IREAD#230) DESC 


~bh 


ORIGSMAL F-A 
OF POOR QU 


KhAC ]KL»Sn#K5lN&#YSI‘^G 

<.!?ITe (i^Sir>3tO) T«iL»5LT# XblNG# YSING 
^gAu (I«£A:),nC) DESC 
WRITfc (UkIT»310) 

DT 3C I = imL#N 

RK AD ( i^EAOf 230 ) XP(I ). YP( I ) 
wPIje (lwRIT»260 ) Xpn)»VP(I) 

3C CQNTlNUc 
L*NL+l 

IP ( YS Yi'^.GT .0. ) GO TO 50 
READ (Ik 6AO>23G) D6SC 
WRITt (IWPIT#320) 

DO AC ] « i» 'JL 

READ (lRtAD»220) VAL»CUM 
► KITE (URIF> 260) VAL»DUM 
J * L - I 
XP( J)« VAL 
YP{J)»i)Uh 
AO CQNTHUE 
GO TO 7J 
50 J»L 

DO 60 I=NL»M 
J*J-1 

XP ( J ) = XP ( I ) 

6C YP u )=-yp( I ) 

70 CQNTIivlUb 

kiPITE {lkPIT» 200 ) Z5(K) 

WRITE (i^RIT.250) TRL# SLT# XSING# YS ING 
WRITE (IwRU>2A0) 

80 CONTINUE 

SC AU = CHG»0 / ( XPd )-XP(NL) > / ZONE 
C(K )»EXP(-PAWZS(K) ) 

XtE(K)»AL/ZQNE + (XSING-XP(NL))YC(K)* THICK* SC ALE 
YLE (K) = YL/ZONEM YS ING-YPI NL )) *C ( K ) *TH ICK* SC A L6 
XX=XP(NL)f(XSING-XP(NL))*C(K)* THICK 
YY*YP(NL ) + ( YSING-YP (?li) )*C (K)* THICK 
CA=CDS (ALPHA) 

SA»SIN (ALPHA ) 

XC.( K) -SCALE* ( aX*CA-J-YY*THICK*SA ) 

YG ( K)»SCALE*( YY*THICK*CA-XX*SA) 

00 90 I-1»N 

XS ( I/K )«SCALE* ( XP ( I)*CA4-THICK*YP( I )*SA) 

YS( I#K)*5CAL£*(THICK*YP(I)*CA-XP( I)*SA) 

90 CONTINUE 

SLOPT(K) «TAN(ATAN( THICK*SLT) -ALPHA )*C(K ) 
TRAIL(K)=ATAN(THICK*TAN(TRL/RAD)*C(K)) 

NP(K)»N 

CHORDO-AMAXl (CHOP DO » CHORD) 

IF (YSYN.LE.O. .0R.ALPHA.N6 .0. ) ISYMO-0 
WRITE (IWRIT>210) ZS(K) 

WRITE (URIT#260) XL » YL » CHORD# THICK » AL 
K«K + 1 


69 


IF (^,Le.NC) GO TO 10 
XK«)(K*((i,-PA) / CHCROO + PA ) 

IF (P4..-<£,1.0) GC TO li?0 
DO Uo l = liN 
CO JCJ 

X(K)»XS( I>K) 

Y{K)»ATAN2{yS( I#K) »ZT(K) ) 
Z(K)«ALGG(YS(I>K)**2^2T(K)#*2)/2. 

100 CONTINoe 

Ei«i. / (Cos( Yd) ) + T1AS iN(r(i ) n 
E2*I./(CLS(Y(NC) ) + T2#S IiMi Y(NC ) ) ) 

F1»(CGS( Y(i) )^T1-SIN(Y(1) ) )+El 
F2-(CJ‘>( Y( 'JC))»T2-SIN( Y(NC) ) )*E2 
El«Sl*Ei 
E?«S2»t2YZT{NC ) 

CALL iPLlF d.NC»Z#)(#A»3,C»i#Ei»l>E2#C»0.#lN0) 
CALL INTPL ( l»NC#ZS»0#l»NC»Z#X»A,a,C»0) 

CALL SPLIF (If NC» Ztr*A*B»C«i*Flf ltF2f C»0.# INO) 
CALL INTPL ( If NCf ZSf Ef IfNCfZfYf Af a#C»0) 

DO llu K=lfNC 
XS( IfK)»D(K) 

YS{If'^)*£IK) 

110 CONTINUE 
120 CONTINUE 

00 13vj K = if.NC 
X(K )=XC (K) 

Y(i< )*ATAN2(YC(K ) fZT(K ) ) 

Z (K)=ALQG( YC (<)^*2+ZT(K)»*2) /2. 

130 CONTINUE 

FA»CDS( Y( 1 ) ) 

EB«SlN(Y(in 

EC«CuS(Y(NC)) 

EO-SIN(Y(NC ) ) 

El-Sl/(cA+Tl*EB) 

E2*SZ*Zr (NC)/{ EC+T2*EC) 
Fl-(EA*Tl-£b)/(EA+Tl^EB) 

F2»(EC*T2-E0)/ (EC+T2+E0) 

call SPLIF ( If NCf ZfXf AfBfOfifElf lfE2fCf0.f INO) 
CALL INTPL (If NCfZSf XCfltNCf ZfXfAfBfDfO) 

CALL SPLIF (IfNCfZfYfAfBfOfifFlflfFZfOfOtflNO) 
CALL INTPL (IfNCfZSfYCfifNCf ZfYf AfBfDfO) 

DC lAO K«lfNC 
X(K)*XLE(K) 

Y(K )»ATAN2(YLE(K ) fZT(K) ) 

Z(K)*ALOG( YL£(K)**2+ZT(K)*»2)/2. 
lAC CONTINUE 

EA»C0S(Y(1)) 

E8*SIN(Yd)) 

EC*C0S(Y(NC) ) 

EO»SIN( Y(NC ) ) 

El-Sl/( EA+cB*Tl) 

E2*S2»ZT(NC)/( EC+ED*T2) 


Of’ 


-7C- ORIGINAL PAGS B 
OF POOR QUALITY 

)/ (t‘A + E‘^^Tl) 
fZ‘(hZ*Tt-E0)/ (ECf£n*T^) 

C4LL SPLIf (l»NC,Z»X>A»3fi3>i.P.l»l^fc‘2,C»vj.#ISC) 

C-iU IMFL a»NC>ZS»XLE»l»M:»2»x,AH»D»C) 

CALL SPLIH ( 1 » N C » Z » Y t A » ri » i) » i » F 1 » 1 » F 2 » D # 0* » I ^ C ) 

CALL INTPL U»NC»ZS/YLF»l#NCf Z» Y»A#3»U»Ul 
G3 TD 1?0 
15C CnNriM'E 
zo«zs( i i 

CO It: K«i»N’c 

Zi (K )«2S {K)-Zy 

160 C Ul^ T 1 N uf; ^ 

17C CTNTlNc't 

cc i';c K«i»NC 

D-J UO 

Xi ( 1,1* ) » A : U »K )-XC ( K) 

Y5( lf< )*YS(I»< )-YC U) 

IPO CONTINUF. 

xy5s^<p(XK♦X5( 1,K) )»CC;S(X»<#y$(i,K))-l. 

Xx=50KT(b>^(2.*XK4XS( l»X) )-^.*XX j-1.) 

XX '>*(XX*XXS)/2, 

xrEo=AVAXl(XTEC»XXS ) 

19C CONTlNud 

ZTIP*Z5(NC) 

200 FDJti'iAr (i6HCPPDFUe AT Z » ,F10.f>/lyH0 T c ANGLE »15 h T£ .>L 

1GP6 »15H X iiNG y-SIMG > 

21C FOPMAf (27HCSECTI0N OFFISiTIDf^ AT Z ■ ,Fi^.6/i6HC XLE pI'j 

IH YLE #1S>H CHORD »1 3HT ^ IC-< N S RATiatlSn al 

2PHA ) 

220 FDRf^AT ( PFiO.o ) 

230 FORMAT (10A8) 

2<fO FORMAT (IHl) 

260 FORMAT ( FIZ .9. 7F15 ,9 ) 

260 FORMAT { faf 15.5) 

270 FORMAT (/5X^HZSYM,12X3HFNC»i0X6HSWFEPl,9X6HS«fcEP2»9X5HStsEEP»10^6«,./ 
lIH£Di»9X6HDIHED2#lCX5HDlH60/lX»8E15.5) 

260 FORMAT ( / 5 X5 HZ S ( K ) , 1 2 X 2HXL » 13X 2HY L > U X 5HC HOP 0 » 1 0 X5 HT h iCR » 1 2 X 2n AL » 1 
12X^HFSEC/1X,7£15.5) 

290 FORMAT ( / 6 X4HY S Y M» 1 1 X 3HFNU » 12 X 3HFNL / 1 X , 9E i 5 . 5 ) 

30 0 FORMAT ( /6X3HT RL f 1 2 X 3HSL T# 11 X5HX S IN3# 10X5HYS I NG/ 1 X , 1 1 5 . 5 ) 

310 FORMAT ( /5X5HXP ( I ) , 10X5HYP ( I) ) 

320 FORMAT ( /6 X 3HV A L> 1 2X 3H0U(1 ) 

END 



Wbiiw"! ••i ‘C - ■» 7f i> ' 4 

Of 


i*Jpl"l"jMNL vv’J^O (NyiNy'*SZ»XTcO»ZT 
iL7iAi»AT#AZ»IY>*^»Ao»Ai#4t»Ai»BC » 0 1»H<£»A, 3»ZfCl »C?»C3» t A » C ) 
bfTS VP i,]i>PTChtb AND j?iN-ii£ C fib? 0 1 N A T(- A 

DIM»\SIJN 4111), 41(1), 4>(1), 43(1), •?0(l)» Pl(l)» y^tl), 
iZ(l), CKl), ..<^(1)1 C4(l), w<»a), C!7(1) 

0 n » ? . / f. K 
C Y * 1 . / ^ Y 
K Y«\YY i 
Z 0 * ^ • u 
K?»N2->c: 

b ; * 1 . / z 

K 1» j 
A H * 1 , 

AY*,!; 

A Y s i . 

AZ-t*- 

*; Z « j . 

»rA X*. 7 i 

IF (H^Y-i.e* J. 1. ) <f^AY«o.h;^b 
S73»S j-^r(7j. ) 

Z'lA X « .o25 

IF (K^yp.tJ.:)) zf-4x»i .-i,/NZ 
SCALsX fLu/ ( . 3J001 »XM4X#xnAX ) 

SC4L2*7iIF/( 1. CCCOOI^Z^IAX ) 

T » 2 . / S C A L 

SY-SUK I ( T ) /TY 

ax« A,>U 1 ( (,(. , ( T-0 ,07 )♦ Aoil ( T-0.&7) ) 

ex=(PtC>/'3-ll, )/F.ft^X)/(F + l.) 
v?» (Dx/DY) 

Wl«l./iCALZ 

WF» ( W1 YD> /OZ i^*2 

bBX«-BXYSO»<T <3 .♦{7.Y-S7i))/((l.+S73)*Xf1AX>*3) 

AAx =l.-&t}X + SQ-i T( (7,YS7i )/12. )*XMAX**3 

CriXs U9,fS 73 ) YX''AX*XMAX/i2. 

AP3X*AHX<-F:iX*{3,YCBX-A.*X'<AX*XdAX)* X^ A X* X 4 X / SOk f ( C 6 X- X M A X * X M A X ) 
l«X»N> + i 
L X * N » / 2 + i 
iO 3C I = 1»'^X 
0D» ( I-LX )*0X 

IF (l.LO.l) 00*-l . + J . /.vX 
IF (I.bw.''-X) DD«i,-i./NA 
b»l . 

IF ( Adi ( UD ) . or . XPAX ) GO TO iC 

A«CFX-(U)*D3 

AS»^bkT(A) 

C * A e X * Ai 4 BdX ♦ ( 3 , *C bx - A . ♦ DD« UD ) ♦ 00* DO 

DO- AB<*UC4dfiX» AS*DOYf 3 

Dl-AS/C 

C2-dtjX*( CBX* (-6.*CBX4 1«^,*D0*DD )--12.*00**A )*0D/( A*C) 

GO TC 2u 

iC IF (L'O.LT.O. ) P--1. 

A- A d-( (OC-3* Xi AX ) 7 ( 1,-XMAX ) ) ♦♦a 


VP 
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C < V \ ) 

L = ( **1 ' ^ * - . . ) V ( i ' - > ) 

[ '' 1 i I • j * ♦ ,, s I < ( < ■ i r ) / ^ 

' i = ^ M ( .' • »-i. > <■ A ‘ ‘ ) 

[ = - u * ♦ Vi M : ^ A >..*>;<( ' . 4^ «' vt u- ) / ( ( .' v- ;, ) / ' ( j, , - v .< , j ,■ 

iMn-.,. 

/*,( I ) = .'<../• • 

M 1 ) = 1 r I 

3 C M n ) - . . V, » V 'ri 
C' .’ ** ■' .J s 
c s ( “ - J i p rr 

4 = , . v,t 

; = i\ ' * A Y 

: * ( » v > ♦ - 1 , ) A { i i ) 

L 1 - ^ C / ( ( A* “ + I j * .’ D 
f : ( j J = A y ♦ j j / c 
f L ( j > = . <■ ., i / ,/ y 
i 3 ( J J = L’ 1 ' i 4 V r 

‘iC • M j) =-uyi> ) (.V. I <1 i + C') I 

- 5 J r ( ^ , » ( ? , + .W* ) > / ( 1 1 . 4 S 7 3 ) + ^ -• A X ^ 3 i 

= i . - V ' r ( ( 7 . * ? 7 - ) / 1 4' . ) ♦ z •I i y 4^ ■•>1 A 

C ■' / = ( i. ^ . 4 ^ y j ) * 7 s A y ♦ 7 j*; t.%f i.i » 

A A A 7 = >■* ■’ 7 4 j } .. M 3 t y" C O 7 - «• • r .' 4 A 4' Z 4i y ) 4'- z ■■• A A “ / * A ‘ ‘ f ( i :« 

•W - « / f 1 

V i. K = £. » A, ^ 

r . ' = ( I' “ IS i ) 1 ) 7 ~ 

1- (K.t,.,.w) •.['»ro-;' . >^:)Z 

(' « i • 

i 1 3 j ( L. »' ) • c r • z i'^ i X ) c j 1 .} j u 

A = C i) 

ASsSCi-TlA) 

Ca»\f;?^Ao4f W 7 ( 3 , ♦ C 3 7 - A , * 0 L’ Y L D ) 4* -J l) * L 7 
r. 7»AcZ’yOL+>i'ZrASTOD<'l‘3 

Ul»zs/C 

C!Z»ijP Z ♦ ( C iZ'!' ( - 6 . dZ + 1 7, 4 JU^uD )-l2. *:j j*4'i ) ♦ VV / ( A*: ) 

GU I - c 3 

( IF ((. U.Ll.j.) b»-i. 

A*i:,-((ub-(j4i7'!4A)/(i,-Z'1AX>)44i3 

c»A4+az 

C * ( A 2 4 A Z-l • ) 4= ( 1 f-A ) 

Cv»« ti4 z iAX4 a5l<z ♦ { ijd-3*z;^aa ) /c 

C] »A*C / ( (1 .4L ) ♦AoflZ ) 

UZa-( A Z+A 7 ) 4= ( Zi'lA A !♦ ( 3 . >D ) / ( ( 1 . 4C ) ♦ A4 ( i . - zv-AX ) ■M'Z ) 
^0 Z ( K J = jCALZ *00 

C L I K ) = . 3 * 0 i. * w i / 0 Z 
CZ ( k )a0i*0i*W? 

C3(K)x. 3*DZ*D2 
C5(K )«t <P(-PA*2 (K) ) 

7 0 CA{K) = C3(k)*C:S?(K) 

7( k2 + i) .Z U2 )4Z (KZ )-Z (k2-l ) 

PFIU^N 

bNO 
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m:..- fy 

QE POOR QUAUW 


<i I^*D# 2U^^ ) 

C GI-.Ntf' <iTt5 >JNGLL4C- LINI- I-J- r« ANSF-C A T1'3N 

Dl-1fcN^joN 2 S( 1 ), XLni)» yLbU)» XC(1)» X/(i)# >22(1 )» fi.(l)» xMl 
1)» iU)f ci(l)p t3(i)» CA(1). bjm 

rn 1 C -< = i,Nic 

K^(K)=«. 

B (X )» :. 

1C CONTI Hut 
Kl«3 
K^ = ^Z + ^' 

KTti=<i 
Cl 2C Ksi^ltKZ 

IP ( 2 ( A) .L r .2V(1) ) KTPl*K*l 
IP (2(A).L5.2S(NC)) KTfc2»A 
20 CONTINUE 

a»CHC(»tlu/2'lNt 
SI « TA^4 ( EP 1 ) 

S2« tan ( S^yt:eP2 ) 

T1«TAN (ClHEDl ) 

T2«TAN(0IHe’J2) 

IF (PA.Nfc.i.O) GO TO 30 
F 1«1. / (CCS ( YLf: (1 M-PTl+SIN(yLE{ in J 
F2«l,/(COS(YLE(NCn + r2»SIN(YLE(NC))) 

T1»(CG4>( YL6(1) )*TI-SIN(YLE (i) n*Fl 
T2»(CDS( YL E( NC n + T2- SIN (YLc (NC ) n^F 2 
S i«Sl>Fl 

S?«'S2*F2<‘EXP(2 S(NG ) ) 

30 CONTINUE 

CALL SPLIf (l»NCiZS»XLE»ui»c2#c3»i#3i#l»S2»C#0»»lND) 

CALL INTPL (Krfl»KTt2»2»>C»^»NC^ZS#XLE>ti#e2» E3.0) 

CALL n rpL (KTE nKTe2»2>X2#if NC>ZS#£1»E2.E3»£A»0) 

CALL INTPL (KrEi>KTE2»2»XZ2#l»NC»ZS»62»E3*6A»t5»0) 
call S^LIF (l»NC#2S»YLE»El#£2#E3>nTi»l»T2»C»C«#lNC.) 

CALL INTPL (KTEl#KTE2»2»YC»i.»NC»ZS» YLtjti»E2» £3»0) 

CALL INTPL (KrEl>KTE2»Z^YZ»i»NC»ZS#ti»c2»E3»cA,C) 

CALL INTPL (KTEi#KTE2»Z»Y2Z#X*NC»ZS»E2»£i»E4*6S»C) 

S»6*TAN( 5‘tvctP ) 

Sl*8^Sl 

S2-a«S2 

T«3^Tan(UIHED) 

T1»b«T1 

r2«e^T2 

XC(2)»i.*(XC(3)-XC(4))fXC(51 
YC(2)«3.*( YC( 3 )-YC(4n +YC{5) 

N»KTE2^-1 



ORIGINAL : i 
-7<.- OF POOR QUALI'i 


le (^.or.K^) GO TO 50 
JO <»C 

iZ- ( 2 (k)-Z (<T 5? ) }/i 
^•e *P(-ZZ) 

> C ( " ) • * C ( •< r f* 2 ) * S • Z Z ♦ ( S 2-i » • t I . - A ) 

v:i^)»TCUT^2)^T*ZZ^IT2-n*U.-Al 
xZ(*»)»(S«-(S2-S)*A)/9 
YZ ( ^ ( I ♦ ( T?-T ) *4 ) /fl 

* 7Z (K » (5 ^-S ) *4/ ( ) 

tZZ(k»«-(T2-T»*a/( fa*d ) 

*•0 CJ^Tl^ut 
tC CONTi.MCt 

YC(^l-n•rC(^l♦l) 

YC ( kZ^I ) ■ YC ( ^2-1 ) 

•i <? T U** < 
k 40 


Sc6‘'C'UTiNfc it'^F (NJ,Nft,NC»NA,NZ»ISYM»KTti»<rc2»SrALi4C»4:#Zi» *C#YC» 
lSL0»*T»TWAlLt)iS»YS»NP»lTtl»ITE?#IV»S0»XP#YP»01»02»03»»»Y,iS0»X'<,,>A, 
2 aZ» YZ» aUCW'^SY'I) 

iNTt-PULATfi iAPPEO blADL SuPFaCE AT **EiH pJInTS 

INTEPP jL AT ICN IS LINEA.< POP CYLIMOPICAL CQOSOlNATtS 

Clf^ENjlJN S0('«E»i)» XS(ND»1)» YS(N0»l)» ZSU)# SLCM(l)» T-AiL(l)» 

i xcm# Yjm» Aom» zui# xm# y(ii» xpcd# ypid# ohd# jzu) 
P, 03(1)» lv(NE>l)» NPdi# ITEKD* ITE2m» xZ(D* fiilif Aid)* C 
3id ) 

P I« 3.1 <.159 2b 5 
Sl-.i*SCAL 
T-l./Sl 
0<«2. /NX 
LX«NX/Z*1 
•ix»\i ♦! 

PZ«N/d 

I VO-l-ISYM-I SY*1-ISYM 
lVl--l-IbYf1 

00 IC K-lf'IZ 

1 TEl (O-MX 
I TE2(K ) -MX 
DO 10 I-l»iX 
lVd#d--2 
SO( I 

10 CONTINUE 
k-kTEI 
K2-1 

20 K2-K2^1 
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OF Q.!AL1»T 


*• l» * ?-i 

i. 

!►' ( ? S(^2 )-Z (< ) ) 2C.<»C.3a 

30 P2-(f»P(^('0)-E*'''l2S(Kim/(6X<»(2S(<2)>-6»»»(23<Kl))J 
ii2«»4*52*(i.--*A)*IZ(*<)-Zj(i<i))/(2i(<2)-Zi(<l)) 

<.0 PI-1.--2 

C«Ci*<S(i.<i)^P2*»b(l»<2) 

0»*i l*r i ( 1# <1 ) S( 1 »<2 ) 

C»iUPT(tJ'f»(:.*XK*C)-2.*C*-l.) 

C • ( C^C X ) /2 . 

CC-3tPr( (C^C)/SC*L) 

03 bL I-2»>lx 

IF ( { AC ( 1 ) ♦ . 5*0X ) .LT.-CC ) 

IF ( ( AC ( 1) -.5*0x ) .LT .CC ) I2«I 
50 CCnTINLc 
ITEi(K ) -II 
ITE2(0«12 
CC»At( 12)/CC 

K K ■ *» 1 
P-P 1 

60 ^^■NP('<^) 

XA-ExP(XK*»5(l,KK))*CCS(XK*rS(l»K<))-l. 

tA«FXP(XK*xS(l»<K))*SlN(XK*TS{i»KK)J 

»0«ExP(<n*x5(N,HK>)*CCS(XK*TS(N»KKn-i. 

Y3-ExP(XK*XS(N»-<K))*SlN(X<AyS(N»KK)) 

C X-Si-^ T ( »A»XA*>YA*YA) 

Cx»(*a^Cx)/2. 

U«StH(CX/C) /CC 
DC 7C I-l.'iX 

X ( I ) >C*AO( I ) 

70 CONTINUE 
AN&L»P I ♦PI 

U* 1 • 

V-0 . 

DO 9C i-i»N 

Xl«ExP{AK*xS(I#KK))*CUS(XK*yS(I»<K))-l, 
yi-EXP{»K*<S(I»KK))*SlN{XK*yS(I»<K)) 
fi-SOfc r ( XI**2^ Y I**2 ) 

IF (P.bJ.O.) GJ TO 80 

ANGL«ANGL^ATAN2((U*YI-Y*Xl)»(U*XltV*YI)) 

U- X I 
V-Y I 

S-SQkr( <P^P)/3CAL ) 

XP( I)-R*CQS( .5PANGU 
YP( I )-?*SIN( .6*ANGU 
GO TO 90 
8C ANGL-PI 
U«- 1 . 

V>0. 

X p ( n ■ j . 

YP( n-o. 


ri 
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wc c:ntini6 

ANGL* AT4N( SL3->T(KK ) » 

ANGL l>ArAN(rA/|A) 

AKGl:-» ATAN(yo/Xi) 

AN(iLi»AS0L-.i*(ANGLl-TeAlL(»k*O) 

ANOL^•Ar^lL-.^♦(ANGL^♦T*«AIL(^*‘)» 

/NGLl«AN(,Ll»*K*r3( 1 »kk) 

ANGL2»ANviL?^*'<A7S(N»*<»‘ ) 
r 1- T AN ( A.nGl 1 ) 

T^•T AN( Af.GLZ) 

CALL $»5Llf- ( l»K»Xr»,YA, Ji, r^,0»0. »1N0) 

call I^TPL ( I 1 * 12* X# y» i» n* Xf>» YP»01 f 32* 03* 0 ) 

D6FIMTICN CF T^^F vQPIFa SHfc6T 

X A ■ > s ( 1 * " »^ ) 

XI • .2^A * A 

A -SLuP T ( ) • ( X A-X 1 ) 
e« I ./ ( X A-X I ) 

anGL I ♦*>! 

L ■ I . 

V » 0 . 

»*5»-l.^jCAL*(x(Il)AA2-Y(m**2)/2,/0/0 
< j AL *x ( I 1 ) AY ( I I ) 

YYS-/rAN2('^YS*XXS)/XK 
Cj lie I-l*.'i 
t » 1 . / 1 . 

ICC caNTiNjt 
t ■ 1 • ( 5 *t 

xx1.aLuG{1.^.5aEaSCALax(I)aa2/C/w)/Xk 
C- &♦ ( X X i-« 1 ) 

YYl »YYS*AA (0-1 . )/0 
XX.eXf>(XKAXXl)ACQS(X-<AYYl)-l, 
YY»exP(«K*XXl)ASIK(XKAYYl) 
p- SwP T < xxa*2a YYAA2 ) 

A\&L»ANGLaATAN2((UAYY-VAXX),(U*XXAV*YY)) 

L« X X 
V- YY 

P-SGkT ( )/$CAL )AvJ 
XP{ 1)-PACCS( .5AANGL) 

IF (XP(i ).GT.x( n ) GO TO 100 
Y( I )-R*SIN( .SAANGL? 

lie c:ni ixofc 
X3- Xi ( Nf KK ) 

A-5LLPr(KKJ*( XB-Xl) 

G- 1 . / ( XtJ-X I ) 

ANGL-0. 

L>» 1 . 

V ■ 0 . 

M-I 2^1 

XXS-1.aSCALA(x(I2» *a2-Y( 12 ) *A2 ) /2 . /0/Q 
YYS-SCALAX(I2)AY(I2) 


77 


Oh . 
OF POr 


YYS »AT 4^2 ( rvs# )l <S I /XK 
C3 130 l-'i.'i* 

£•1 ./i .09 
120 CONTlNOt 
£•1 .09«e 

»<l«4LJG(l.».3*F*SC4L*xn)4*2/C/0»/XM 

c« ( * <i-x I ) 

Yyi-Yrs*A*( 0 -i. )/o 

YY«6>P(*n*«<i»*biK(«K*rr4) 

»«S54T(oy*2»YY*»2) 

ASGL»AN0L4ATAN2((L*Yr-V*««), <U<<X4V4YYM 
U«x X 
V • Y Y 

P -so** T (((.♦•')/ sc AL ) <0 

xP( .94AK0L) 

IP MP( 1 ) .L T . X ( I M GO TO 120 
Y( I ) •><*5IN( .5< ASGL i 
13C CJN I iNtJt 
C-P<CC 

Cl 14C I-I.MX 

S0( I »•<) -see I »< » ♦O* Y( I ) 

KC CI M lNu6 

IP (<K.tc.<2I GC TO 150 

K < ■ K 2 
P« R2 

GO TO tj 

15C Cl 160 I-il» 12 
I V( I »< ) »2 
160 CONTiNUs 

SO (!»•<)• S0{2»<MS0(2#K)-S0(i»<) 
SO(«>>K)-SO<NX,K)4SO(NX»K)-iO(NX-l,<) 

M- I 1-1 

00 170 I-2»i 
ZZ- Z (<) 

IP (ZZ.G£.Z(kTE 1M lV(I»K)«iVO 
170 CONTI'luE 
6 ■ I 2 ♦ 1 

00 160 1-H,nx 
ZZ«Z(K) 

IP (ZZ.Gfc.Z(KTEl ) ) IV(IfK)«lVO 

160 continue 

K2-K2-1 
K» K ♦ 1 

IF (K.Lt.KTE2» GO to 20 
Kl-3 
K 2 ■ N 2 <2 

190 DO 200 I-2»NX 
ZZ« Z (< ) 

IP (2Z.lE.ZS(NC ).AN0.ZZ.GE.2(KTE1M W(I*K)«IvO 
200 CONTINUE 
K»K ♦! 


ORIGINAL P/ 

OF POOR QUAL.i'V 


- 7 h- 


IF Oli Tl» 

N • T ( ^ 
i • I T f 1 ( ^ T r n 
OL i ^ J <•" L > ^ 1 
on :i : i 

IF ( 1 »■( I »a) .or .J) G*J 

?IC CINTI iLL 
b • SC ( L * » " ) 

IF (;,/•»«• J/ ( r .L I . 1 . fc-w3 ) IV(LI»K)«J 
2?C CJ.'.HNIE 

I t 3 0 1 ■ 1 » '’ * 

iC(i»'i/)«3.*(3(/(i.'^:-n-$c(i»f<z-2n*so(i.sz) 

J>3( I»<1-1 ) 'SC( I»^l♦l t 
23 0 CDMlNdc 

ITFl(<l-l)»ITlf KKl^n 
ITt2( iZ ) T6 I (.iZ ) 

ITt2(< i-1) -ITiZIKl^n 

►. f TO' N 

LNC 


bOe^CLT i%t tali'* 

initial tsri**ATt OF VfD'JCkO POTENTIAL 

C1'*'NCn &(12N,2t»17)»SC(U'V»i7)»E0(l7)»IV(l2^»i7),ITt3(17)»ITF2(l7) 
l,AOn^n^ALn^<f) 0 A 2 (L^)),A 3 ^ 12 i),^C( 26 ^,U{?t),t> 2 Uof^•i}^ 2 t))^In 7 ^ 
2#Cl(17)#C2(l7),C3(i7),C<.(17)»C5(l7)»XC(17),*Z(l7)»>2Z(l7)»YCU7),> 
3Z(17)»TZZ(17),KSrr•,^^^,^T»N2,■<TEl»KTE2»IbY^'.SCAL»bCALZ»»<»a'*Ef•A,ALF 
<»HA»Ca» 5a »F » T »PA» At*cAJc 
T«2 , / 5C al 
K r»NY^ I 
FZ-NZO 
on 1C 1-1»129 
UO 10 J- If 2& 

00 10 K-1,17 

1C CONTINUE 
••Z-3 

DO 30 I-2fNY 
00 20 K-kZ,mz 

IF ( I Y ( i ) .L T .2 ) GO TO 2D 
0SI-iv)( I + 1»M-SC( 1-1»K ) 

OSK-Sj (!#<♦! )-S0( I»K-1 ) 



-7^- ORr 

OF ‘ 

* » • A I ( I ) • C j I 
57*C1 

I • a;< I) 
y • b V <1 . ^ ) 

»( (T*«**i-T**i )••?♦(?, 

y* 

H» I , / r ►• 

0Za-ii«<4(<)*r^«iZ(< » 

A -C t ( < ) 

?«A*C 

FZZ- 4* < 

F VZ-H*F^/* ii 
F»Z»»-*FZ4.«a2 

^Yfa^<«x«(riaritFZZA((^Z*dZ^'(A«xji)) 

Fil>M«H«(^i«li4F7Z*(A^«AZ«yxAyX)) 

AiarkT*S>^-<''>-5ZA<=xZ 

4 Ya? YV-i,ap « >-b7*FyZ 

47«Ft Z-i»*F* Z-SZ*FZZ 
‘•yBAY-i>A*A»-5^*AZ 

Y»i4*«A**'/SA7-CA*yx-J'^tv»A*IX/C<»(K) 

U»(,A*»<aS4>y>»P/FZZ-'J^cGA*TA/CA(<) 

• •CA•*^(")♦5&•(C♦P•yZ('<))/^ZZ-01E&4•yZ(^)/C•»(^) 

G( if-'ytifK )«v,( l»Ky-i»<^)Y(4x*!j^Ay*v^AZ*w)/((iy*»-i(Ky)) 
? C w n < T I 'Y 1 1 
30 G ^< M '■'i L c 

K i a n T £ 1 
•< 2 • Tc 2 
03 AC 

e ( ( » ) • 0 . 

AC CO'*Ti\ut 
Tl*' 4 
FnD 


iU-JFOOTlNF ^JjNO 

CFFiNfcS THc eJUSOA^y VALUEb rjF THf Vt'LJCITy PJTFMiAL 3 
C3''*'0s 0(12 J»Zt*17),SC(l?'9,i7), £0(17), lv(12 5. 17), ITEKlTliIT-Zd?) 
l,AD(l?t),Al(i2-i),A2(U^),A3U29),(‘0(2t>),^l(2o),t»?(2o),d3(£o),Z(l7) 
2»Cl(17),CZ(17),C3(l7)#CA(17),C5(17),XC(i7),<Z(17)»XZZ(17),yCa7),y 
3Z(l7),yZ2(4.7),KSyM,Nx»U#N2#KTEl»Kr£2»I3yM,SCAL»SC&LZ»A^#3.'l6i,A,ALP 
<,hA,CA#jA,FiAC*ifT ,PA,A6LA3E 


00 - 


:;,c: 

OF POOR QUALl'n 


/CAL/ Pl»P2#P3*dETA»f>R»BNfOG»GM»NS*Ji*vlfkiI»l*>‘fJ-<><<»lOtJJ» 

IKG 

COHnON /SWP/ Gl( 129*26)* G2U?9f 26 )*SMi2W)»S2(WO)»Sxk(12<}),S(7( 12 
19)»bZ2(129)*ri0(129)*Pl(129)*C(129)*3(129),li.i2*Lx»MX,<T*-r>ri.AAC 
2*31*C2»N« 

L *«N* /2ti 

f<««NX*l 

K Y«NY^ 1 

PY«KY42 

^Z*N2»3 

0X*2«/NX 

U«0x •ox 

AA0«1 . /Ff*ACH**2^.2 
ZSYM-1-k>y« 

C1-2./P1 
02*1. /?2 
PM*0 . 

GM-0. 

NM«0 

Fa-c. 

IB-0 
J R«0 
K Pa 0 
OG-0. 

IG-G 
J G ■ 0 
KG»0 
NS-C 
11-2 
I2-NX 

DO 10 J-1*1Y 
00 10 I-l*'iX 
Gl( 1*J)-G( 1* J*l) 

G2( I*J )>G( 1* J* 1 ) 

10 CONTINUE 
JT« 1 

20 CONTINUE 

KL« (NZ^6-JT*(NZ^1 ) )/2 

KMaKL^JT 

KNaKMt JT 

I«LX 

G( l*l*KU-G( I*1*KN)^ZSYNP(G( 1*2*KL )-G ( I*2*KN ) ) 

0SI*S0( I^1*KH)-S0( I-1*KM) 

0SK«S0(I*KN)-S0( I*KL) 

SX(I)*A1(I )P0SI 
SZ( 1)«C1(KN)P0SK 
P-1.0 

00 30 J-2*KY 

G(1*J*KU-G(1* J*KN)4ZSY)1P(G(2* J*KU-G( 2*J*KN) ) 

G(NX*4*KL)-G(HX* J*KN)4ZSYN«(G(NX*J*KL)-G(NX* J*KN) ) 
YP>B0(J)^S0(I*KN) 

IF (J.EO.KY) P>AMINO(l*IVn*KNn 


I 


OF t 


» ' 




/. » i * - 1 ( : j < j T 
• < ‘ • ( • < ♦ • ( ! ) ) - j I • ' i ( I 

. »* ( 

' )(*»«)*» 

' » I , V . i i « ( : . J ) ♦ : ^ > . * ( ( ; 

1 

w • • < 

w « • » < * 

( 1 » V » ^ L ) * > ( i » J I • 

; ( 1 1 V * . ) * '■ M : , « ) w . f ' » ( . ( I 
• * >•/'-» 

- J J ' i 4 * 1 • ‘ 

J • L » - I 

0 I _ ^ 

L . M ■< I r 

1 * L * ♦ ; I 

I • . r J 

M » ' . ( i ♦ j . ' ' ) - 5 4 . ( i - ; . ) 

f ‘ 1 »•< L ) 

.• • ( J ) * • M : I » ' ^ I 
M 1 ) » C I t ‘ ^ 

*( I •]•<». )« >( ) 

. ' ' 4 w « f . ^ Y 

‘ o I I ) 

Y * - « ( J ) ♦ ' • ( I • ' ^ ) 

I It,, / t-Kirf *( 

■•» 1 ,/(>***<^YX*Y*) 

/ 7 t - 1 < * I / ( < ' ) - y » ♦ y / (< ,1 ) 

.‘li 

J i ( I ) 

* '. * i - i r 1 * ( > ) 

I 1 • ‘ ( i » V t ^ ‘ ) - » .i ( 1 » J # >» N ) 

1 • J*' 

:m i.j)».,( i.j»KN)o.*(j( Iij»kl)-g(i»j»-<w) ) 

V ( 1 » J . : ) * (m' ( I # J ) ♦ Z > r .Y ♦ ( i ( 1 1 J , 1 ) -Gii < I . J ) ) 

C' ^ T 1 S.r 

J » K y ♦ I 

o(i»J»-L)«(Cl(KD*b(I#J»KN)^A*(G(iP#J»'<U^UGl»-t3*0GJ)/(:i(<’1)^A) 
ci(i»j)»o( :*j»Ku 


•w ) 

2 I 

) » • M 1 ( J ) • J T 
) 

) 

-U* L'GJ ) /C I ) 


Iij)) 



ORIGINAL PAu- 

-M2- OF POOR QUALilV 


1 ► ( I . l I . L < I v'* L 1 ' i ^ C 

7 C COM iNjt 

IP ( J 1 .fcw.-i ) GO Tu liC 

M M a 2 

L - 3 
I ■ 5 
j »Nr 
K aL 

X > 6C( i ) 

rC<<)*AT4u/(2,*i»y, T*A*x-r*Y )/>k) 
CPN«(WA*»-r*Y)**2f(2.***y)**£ 

T4-a*x^y*y >/3cS 
Y«a^,/(M»Y«( T-X*«-Y»Y )/0cN 
H»l./( 

P-Ci (-» )*C3S ( Yl ) 

C - : i ( ^ > • j i >1 ( Y l ) 

AZa-«**l^(<)-YI*YZ(t') 

aZ»-X»*YZ(<)^f»*«Z(K) 

S*(I)-Ai(l)*<jC(I^i»K)-i0(l-l»K)) 

S/ID'CUM’ISOd.K^D-SOIUK-l)) 

&Yaf*l(J)*(j(I,j-l,K)-r,(I,J^i,K)) 

GZ»CU'^)«(G(I»J#*'>l)-&(l»J#K-in 
0»G*-5 X ( 1 > *G Y 
V “ G Y 

■ •GZ-SZ ( : ) *0 r 
LO«H*( X x*U-Y t * V ) 

VC»H*((P»YX4j*AZ)»U^(P*XX^3»bZ)*V)^C*M 

Ul" i*Uu 

l*V'J 

OP lOu K «K < , MZ 

IS«ISIGW(1»^?*(«Z-K)-1) 

IT«1^ lGN(i,^*( K-KK )-i ) 

CO 9C 1«Z»NX 
X»AC( I ) 

Y-SO( I»K)4-dO(Z ) 

0EN-(T+X*X-Y*Y»**2*(2.*X*Y)**2 

XX«2./AK*x*(TtX*X^r*Y)/06N 

YX«2./XK*Y*:(T-X*>-Y*Y)/D6N 

YI«PA*(YC(K)+ATAN2(2.*X*Y,T+X*X-Y*Y)/X<) 

G2( I X 1 ) >G1 ( I .1 ) 

Gl( l>i)-G( IflxK) 

G(I»l»K)-G(I»3#K)-(YX*Ul-XX*(Vl*COS(Yl)-Wl*SIN(Yl))/Cb(M)/ol(2) 

G(l»l.K)«ZSYM*G(I»lxK) 

iX(l)»Al(I)*(SO(I^lfK)-SO(l-lxK)) 

SZ( I)-Cl(iO*(SO( I»K+IS)-SO(I#<-IT) ) 
DSlI»3UI^l»K)-SC(I#K)-30(IxK)^S0(I-l»K)^A3(I)*nSI 
DSKK.S0(I»K^IS)-S0(I»K)-S0(i»K)^SC(I*K-IT)+C3(K)*DS< 
csiK-so(i+i»K^is)-so( i-i#<^is)-son^i#K-m^so( i-i< K-iT) 



i 3 - ORirr'^.. 

OF Pi 


SZ2( I j-CitOOSxK 

eo T INOC 

'J ) <ic 

X » A V. ( ^ ) 

! 0 ( j * < ) «• 1 'M J ) 

r»x*x-r*v)/x<) 

C;r\«( 

YI«?,/*^*Y*(T-X*X-Y*Y)/OcN 
P«Ci(^ )*C0S ( Yl) 

G( l»w •*“ )"?5YM*u( l»J»'< ) 

M ■ ^ X - i 

X« 4 v. ( NX ) 

Y■SC•(^^X»^)♦bO(J) 

ri»FA*(rC(<M4TAK2(2.*x*Y,T4X*K-Y#Y)/X<) 

LE\»(T^x*x-y*y)**2^(2.***y)**2 

X«»c,/X<*X*(r^x*x^Y*Y)/i}tN 

yi.p./iKYY*>(l-X*X-Y*Y)/JES 

►'-cicx )*crs( Yl ) 

J»:5(K)*iiY(Yl) 

GY«ollJ)»JG(K*»J-l»K)-G(tMX,J^l,*<)) 

G2 ( ‘'X ♦ J ) »Gl ( MX , J ) 

GKMx . J ) »&( MX f J .K ) 

G(MX,J.K)-G(M»J»K)^(XX*.Jl^YX/C<»('<)*(f'*Vl-G*hl)tSY(NX)*oY)/Al(NX) 
G(MX» J.O-ZSYM^GCMXjJiK ) 

<yc continue 

G{l»l»<)«G(l»2»'<)*G(2»l.K)-b(2»2»K)*ZSYM, 

G('1x,i,<)«G('‘x,2#K)^G(YX,l,K>-G(Nx,2,K)*Z5YM 

t-G(lTE2('‘)»KY»K)-G(ITfcl{'<)»KY,-<J 

GdiKrylj-cj-GCMx.KY-l^'O-E^ZSYM 

GCMX.^Y-xlf-O-Gd.^Y-l.Kj^E^ZSYM 

G(MX,<YfK)«G(l»'<Y»K)4E»ZSYM 

IF dS^lT.GT.O) CALL YS-EEP (<) 

100 continue 

IF (KiYM.NE.O) GC TO 110 
J T»- 1 
GO TO 20 
11 G CONTINUE 

M R « 1 . 2 ♦ F K / A A 0 

R**»1.2*<<M/nm/AA0 

GM«GMZn'1 

RETURN 

END 





ORIGINAL PAL. T 
OF POOR QUALnV 


S'Jb^aL•Il^fc YS4 EcP (k) 

c Th 6 EJOATIJNS f-0« & APE SGLyEU nfePE »F0K Ml«Eu SUdSOMC 

L AnjO SL'PE-<SJNI: PL3* »EY eDn PPLA<aTIJn ,AND By USING A 

C shtaTEL OI-PtHtNCfc SChc'-E 

o(i2^.£t»i7),SC(i<J9»i7).E0(i7),IV(12<;»l7)»ITE4(i7),ITE:(17) 

2»Cl(i7)»C2(i7),C3(l7>,CA(17)fC5(l7)»xC(i7),W(17),»ZZ(l7),YC(i7).Y 
3Z(17)fYZ2(l7)iKSY^,N<,MY,SZ»^TEi»KrE2»ISYr^»SCAL#SCAL2»XH,orEGt,ALP 
<iHAtCAfiA»FMACH»T #PA»ABLA06 

CO^i^CN /CAL/ Pl»P2»P3»8tT4»FS»Fi»OG»G'^»NS»Ul»Vl»bl»I».»j»i»-<fc»iG»jr,, 

IKG 

CC.l'ij.'* /b^^/ Gl(129»2t»»G?(12?.2b)»S<(l2'J)»S2(129)»SXA{i2v),SxZ(i2 
lv)»iZZ(lc'#>»^0(l2^)»«l{l2’»)»C(129)»0(l29)»li»I2»L<»'i»#KY,M/,Ti,AtC 
2 » G 1 » c i , s ^ 

L-H 
J !■ 2 

IF (Ff^ACh.GE . I . ) Jl«3 
C ( I 1-1 )-C. 

D( I 1-1) «C. 

DO 1C i- II. 12 
wj( I )-l. 

Pl( 1 )>1. 

Gl( I.1)«G( I.lfL) 

Gl(l»Ji-l)-o(I#Jl-i»L) 

1C CGNTINUc 
J- J1 
13-12 

2C COnTINJc 

ec--Ti*8i( J ) •: i(K ) 

D3 to I-I1.I3 
AB--T1*A1 ( I ) *81 ( J ) 

AC-T1*A1 (I )*C1 (K ) 

X>AC( 1 ) 

Y»SO( I .K )^.30(J ) 

0EN-(T^X*x-Y*Y)**2^(2.*X*Y)**2 
XX«2./XK»X*(T*X*X^Y*Y)/0EN 
Yx.2./xk*y*(T-XAx-Y*Y)/DEN 
FH-I.-xOCI )^XX*XX4-YX*YX 
H-PO( I ) / FH 

AZ--Xx*A2C<)-rx*YZ{K) 

8 Z • - X X • Y Z ( <) ♦ Y X ♦ X z ( K ) 

E -H*AZ 
F«H*8Z 

XXX«2./XK»(( T^3.*X*X^ Y*Y )*OtN-A.* X* X*(T^X»X^Y*Y)**2)/OEN**2 
YYY-A,/XK»X»Y*(0EN*2.*(T^X*X^Y*Y)*(T-X*X-Y*Y))/DEN**2 
AA-PO( I ) «C S(H ) 

Y1-YC(K)4-ATAN2(2.AX*Y»T4>X*X-YAY)/XK 


•4 


-► 


= Jh , 

OF Tl ^ 


' I 


V. I* T‘‘C rwLL.wl^*^ LiN-S *r-*t Cu^^^lCIkf. j*-' T>-t 

W •* f J • i I »Ai «hLL AS Th*I^ 

c \>. th-: c JL-": r. A I c s 

£ ■ H • ( I 
>•■►«• T I 

)*YVV| 

^. • A A « j i V I Y I • u ) 

H » « - . • Y t 
><<•— .»! ^ 

PY»— ► 

••(■ 

t»»-(i»^»;(* )*AY•r7(^)) 

M»-(AT*«/(')-A»*Y^(-<)) 

h7»-(A**iJ(-<)^c*Y/Z(»v)> 

P£«-A*Y4/(<)*3*|//(K ) 

YYJ»r*A^.'*P 

y t ^ mr i -4 -s* .Y"t->>*AX4';*6Y 

YYV«PY*A^'. r*F<-P*AYti*t'< 

YY«»H7*Afw7vFfC*FZ 



* W 

< <f 0 ♦ 

P*E*- 

U*4 Y 













f xv*>'Y*e 


Y* 

P •£ Y-f 


A X 













/»•.>» ^ Z’P 

• W 

7 * A ♦ 

^♦tz 















7 Y ^ • > • A 

•* 4t 

< • £- 

p*e Y- 

c* 

A X 













i f \. - P Y • P 

» 

• 4 

Y* Af 

p♦ c X- 

«:♦ 

AY 













tY,<«PZ»r 

• W 

1* U* 

p* P z 















CGl»o( !♦ 

if 

J f L ) 

-G( I- 

1 • 

J f L ) 













Cr.j «r ( I , 

JY 

I f L ) 

-G1 ( I 

f J 

-i ) 













C G ^ ■ o ( If 

J f 

] ) 

-Gi ( I 

f J 

) 













roi ( 1 

♦1 

f J f L 

)-G(I 

f J 

f L )- 

G(1 

f J 

f L ) ♦ G ( I 

- 

1 

f J 

f L ) 

Y A3( 

I ) 

* 

)v.’ 1 


CGJJ“o( i 

f J 

♦ i f L )-G ( I 

f J 

f L )- 

i( 1 

f J 

f L ) «>G ( 1 

f 

J 

-i 

.L ) 

-c3( J ) 

• 

L) w J 


UG^ < « G ( I 

f J 

f L^l 

)-G( I 

f J 

f L 1- 

G( I 

f J 

f L )Y-G( I 

f 

J 

f L 

-i ) ♦CBC 

K ) 

• 

0 G f 


L 0 1 J * G ( 1 

♦i 

f J ♦ i 

f L ) -G 

( I 

” i f J 

♦ if 

L ) 

“O ( I ♦ 1 f 

J 

- 

if 

L ) ♦ 

Gl 1- 

i f 

J 

” i. 1 

L ) 

CGI K-G( 1 

♦1 

• J f L 

♦ 1 ) “ VJ 

( I 

i f J 

f L- 

1 ) 

-j( I-if 

J 

f 

L«-i )♦ 

G ( I- 

if 

J 

f L- 

1 ) 

JG j »G ( i 

f J 

♦• 1 f L 

♦ 1 )-G 

( I 

f J-1 

.L^l ) 

“G ( 1 f J ♦ 

1 

f 

L- 

i )♦ 

G ( I f 

J- 

i 

f L- 

i) 

GX» £ i ( 1 ) 

• u 

J 1 
















GY»-rtl ( J 

)♦ 

J J 
















G2-C 1 (^ ) ♦O 

3k 

















X ( I ) *&Y 

V«-GY 

•y-'-GZ-SZ ( I )*GY 

U«U!i^:A*XY*SA*YXAO/C<»(^) 

V«Ve^SA*X<*P/CA{K)-CAAYX 





ORIGINAL PaGu ‘.3 
OF POOR QUALITY 

- j* ^ 

v0« 

I i J*J* 7 1 : * ^ **'* m 

V “)* VO-u"t jA*PKh( »<) 

F •j^YTJ*Trj^?fJ*Z>J 

rYY«-a«^tT»J*YYJ^Z*J*/YJ 

F»/«', 

FY/ «l.*TYJ4t»*ZYj 

iV«-b*U^^YYj*VO^ZYJ*,& 

P*«0*YC+F*WL 

CU»A*AX-t:*AY^YXJ*YXIJfYfJ*YXwt3*VX*^ZxJ*ZXU«’7YJ*ZXVt-*i*.. 

0V»-i*4Y-fc*iX^V«J-»YYv^rrj*fTV*.*YY,4Z*J*Z»J»ZYJ*ZYV»F*Zt# 

0-»r>J*^X4-YYJ*CY4':*wZ^Z*J*PX^ZYj*PY4>P*PZ 

Cl»bO*(AX*>'JtAY*»iV)^VC*(Y*0»PU^YXV*^V^YX-x*ti-.)4’ixJ*(ZXL'*d'."*ZXV*;?V^/« 

1 i* • 3 K ) 

CV'wO^C-AT^^C^AXYtJYMYj^lYf^^aU+YYv^’jY^YYn^oKl^wC^IZYL^ol^ZfV^cV*/ 
i Y ** r w ) 

C-»wC*( - *•! 

F X X - H . ♦ « 1 ( ; » - . 0 ( 1 ) 
f Y f .F T r ♦ 1 .-- ; ( 1 ) 

FYY«S<( 1 )**?*rXX4FYY«-SZ( n**2*F 
1 I )♦S4( i )*F X^ 

F*Y«Fxf-i,Z( I)*Fxa-3*( I)*Mx 
FYZ« f YZ-i < ( I ) •FXZ-SZC 1 )*PZZ 
Av»fcv-b * ( I ) »ij-sz ( 1 

UU*(Ju*rXo 

I Y « •: t • A V 

U ^ U ♦ " 

VV«AV*AV 

V 4 • AV ♦ 0^ 
i»W»t?l4*e<X 
AXX«FXX*A A-li'J 
AZZ*FzZ*AA—»« 

AXZ«a. * ( Fx ?♦ A A-UW ) 

F»-UXA»iX<(I»4AZZ*SZZ(l)4-AXZ*5XZ(I))*Gr*n*l ( AA • Jl-C c) •tt‘4(AA*.!V- 

IC Y ) •Vti^ ( AA YOn-C *» ) ) 

AXT«£3S(bL'YAl(l)) 

AYI«A3S(AV*Ml(J)) 

AZr«AJS(b4*Ci(tx) ) 

A«-xj(i)*rt:rA»AA/A»-AXi{Axr»AYT»AZT»l.-fc)(I)) 

axT«a*axT 

A YT » A* A Y f 
AZT » A» AZT 


X«L'U^PY*jY + PiYbW) 

ZZ-Z .• j X ( I ) *F X >-2 . * 'Z ( I )• F rz ^ ‘ ( 




JJ - OR!c.ii . I .4 

OF POOR C. / 

oU’<>3nIC •*JI'»Tb Sc'^AJATtL 

(«a»OC«'^<^) C'- iC 
£k»aAi(<»4^(I) 

A r Y ■ ( i: y T i- V V ) *b2 ( J ) 

A7. 

A*Y»(*:if*Ai-uV )*(A^tAf;) 

A*Z*A»4*AC 

AY7«(rt'*AA-\< - » 

c I* ■ A > » 

•jM* A * * 

< jyy^aZZ ) 

^»■A»^♦OL^l♦AYY*OOJJ♦AZZ♦CG•<^♦AKY♦^^lIJ♦AY^•CbJ*<.♦AX7•OOI•^♦^< 

GO TC <*0 
3C Cr.NTINcc 
N S ■ N j ♦ I 

]M. 1- iF i * ( J 1 ) 

!'***« M-lFl X ( 31 ) 

A * » «UL»* A z ( i, ) 

AVYm^>/#e2( J) 

aZZ**>**C2('X) 

AXY.‘j.«>i*J\»*AF 
AXZ«^.*S1*J.«*AC 
A YZ • V •• 3C 

3t«a(F\X«JwOJ)*Ae(I) 

FYY»(FtY*UJ-i«Y 
bZZMFZZ* J J-Nd )*C2 (X ) 

B«y»(Fxt*^j-UV)*(AB^ah) 
bXZ- ( F xz *G w-L- ) ♦ ( AC^AC ) 
bYZ«(FYZ*'JJ-V'a )*(bC^bC) 

AQ>AA/Gw 

GkLrAG»i**YD3IlYBYY*t;GJJ^3ZZ*0GXXfexY*0GTJ*6YZ*DGJK4.t>*ZYGGl< 

DGII«G(1>J»L)-G{I.'1,J,L)-G(I.'1»J»L)+C(IM'1,J#L)ya3(1)*DGI 

LGJJ-G(I»J»L)-Gn#J-i»L)-G(l»J-l»L)^bl(l.J-Z)-P3(J)*t;GJ 

DGKi<«i,(i»J#L)-G(l»J»L-l)-G(i»J.»L-i)*GZ:(I»J)>C3(X)*0GX 

CGlJ«G(I#J.L)-G(IM,J,L)-G(I»J-l»L)*G(Ii»J-l»L) 

CGI<»G(i»J»L)-G(I»J»L-l)-G(iM»J»U^G(I'1,J,L-l) 

DGJK.G(i»J.L)-GU*J»L-i)-G(i»J-l»U>G(I»J-i»L-l) 

GS3«AXJi*0C-IlYAYY*0GJJ^AZZ*0GK<^AXY*3GlJtAYZ*0oJK*AXZ*CGlK 

=«.5*(AG-i. )*( axx^axx^axy^aaZ) 

pp« AC*bx *-( 1 .-S I) *d 

bM» A3*3 X X- ( 1 . 4- S I ) *0 

C,«-Aa^(tiXX4>bXX4a2Y(0YYt3ZZ))4(AQ-l,)*(2,*(AXX*AYY4AZZ)4AXY4AYi4axZ 
1 ) 

w»(AO-l.)*GSS4-Aa*0ELTAG4P 
AO CONTiNOt 

IF (AdS(«) .L6. A9S(FP) ) GJ TO 50 

FP>k 

IP. I 

Jp ■ J 

KP »K 

50 CCNT1NU6 


r2 



ORIGINAL FAu-; ‘S 
OF POOR OUALirV 


-8f 


N M ■ \ ,i ♦ i 

3»d -4*l-AfT-A/T 
e^** fa A « T 

fa ■ I . / ( I i- I ) » 

Cl 1 )»3* iP 

&l I J-3* t ►-8^*0 ( I-l ) ) 
bC CTNTiNLt 
C G ■ u • 

1-13 

00 ec r.ll.13 
cc-LM 1 )-ci n*co 
G^'■Gr'♦• A-j i ICG ) 

IF ( A3 i ( cG I . L-: , AfaS (OG ) ) GO 10 70 
w G ■ C G 
IG- I 
jG> j 
K G ■ 

7C CJNU>IU£ 

G2( I«Ji>Gl( 1/ J ) 

GI( 1>J )>G( 1» J» U 
G( I#J»L)*Gll»J»L)“CG 

1 -I-l 

fO CnKTHct 
J» J* I 

IF (J-KY) £J»-?C»110 
90 CONllNtt 

IF (I3.GT.ITfc2(«)) I3«ITt2(K) 

IP ( ITe2(t<) .Fw.f'Ji) I3-LX 
CO 100 l«Ii.I3 
LV- Udfa ( 1-IA3S I I V ( I»< ) ) ) 
kO( I)>AnINO(LYtUBS(IV(I»K))) 
kl { I ) -LV 
IOC CONTl'JUE 
GO TO 20 
lie CONTINUE 
I-L X*i 
IO-Nx^2-I3 

IF (K.Gr.KT62) GC TO 130 
DO 120 1-I0» I 3 
X- A0( I) 

Y-SOI ) 

0EN-(T*x*X-r*Y)»*2+(2.*X*Y)**2 
XX-2./X-<*X*(T^X*X^Y*n/0EN 
YX»2,/Xrt.»Y*IT-X*X-Y*Y»/0eN 
A-1 .-kO( 1 ) >X X ♦XX + YX^YX 
H-SfC(I )/A 

A7»-XX*X 2(<)-Y X*YZ IK ) 

BZ.-X X *YZI K) ♦Y X*XZ IK ) 

A»ftOI I)*C5IK ) 

Yl-YClK)4ATAN2C2t*X*Y»T^X*X-Y*Y)/XK 


^ 7 Z • A* A 

Fvr.H*^#(Y<*Y»^PZZ*(^Z*'»Z*«AAJi*))yl.-’<0(n 
1 1 1 • ( I <•*! ♦^zz*(AZ*AZyyA*Y«) )y4»-^C( I ) 

FlYaH*i4«(>/k«YA»FZZ*(A/* W^aX^Ya) ) 
f t l» i ^ > \? 

f-f il* U 

AK«F»Y-j*( I)* r><-xZ(I)*P*? 

AY»rYY-i»(i)*CAY-S7(I)*FYZ 

AZ'FfA-xCDyj'AZ-sztn^'-zz 

*iY«AY-:A(l )*AA-jZ<n*A/ 

oGi»c(i^i.<Y,L)-C(i-i»'^y.u 
'Gi I • »»- Y , L ♦! ) -G? (lp>^y) 

i-CAyyx-j*t*, a**a/C<*( - ) 

c-Al(I»*uoly:A*AA^S«*y**‘"/G<*('<)-G'’EGA*YA/C<*('‘) 

• «Ci(<)^L3''^CA*»Z(*‘)ySA*0^»>*YZ(<)J/C<*{<)-G''PoA*YZ(<)/;<.(<) 
G(I»^Y^;.L»»i(l»KV-l#L)^(AA**J*AY*VyAZ*w)/(rY*-l('<Yn 
123 CG'NTiNUt 
I«IC 

IP ( i:.Nc . I TFi (< ) ) GO T3 1 3w 
E-Stnt^YfLJ-vjCIOf^YfL) 
c’ClU-:J(L)**»3*(t -60( L ) ) 

130 CCM iNuE 

IP ( I . Lc . I I ) i F TUKN 
I -I-l 

c ■ w • 

IP (1v(4*‘*)*ne*i) go TC I've 

t »E 0 { L ) 

1<*C CJNTi^Ut 
»'«na^ 2 -I 


G { I, 

KY ♦ 

i A 

L ) • 

G( >*t 

K Y-l » 

L ) 

-6 

0 ( ^ A lA Y ♦ 

i A 

L ) ■ 

SJ 

( I A 

X Y-l, 

L ) 

♦ t 

G2( ^a-a Y 

) « 

Gl ( 


A K Y ) 




ol( 'I 

aiA Y 

1 • 

G(« 

A lA Y A 

L ) 




G( Ma 

K f . 

L ) 

«G( 

1 

A *« Y 

A L i 





G3 TO 130 
ENO 


S'JBPJoTINL vELO (K»L*SV»5^»CP»a*Y) 

CALCULATtS SU-'FACt VELjCITY aNC P^ESSL'^E COlPFICIEM 
C0'1^':N o(i2?»2t»17),SC(i2^#i7),fcO(17),IV(12«ifi7)»ITti(17)»irPi(17) 
l»A0(i2Y)#Ai(i2'i)»A2(12<;)»A3(12'A)»BC(26)»'<i(2o)»32(2o).'J3(2t)»Zll7» 
2.Cl(17)»C2a7),C3(17),C‘»(l7),C5(l7)»ACU7),xZ(17)»AZZU7).rC(i7),Y 
3Z(17)»YZA(I7)»KSYM»NX»fyY»NZ#'<T£l»KTP2fl5>Yi1»'5CAL»SCAL7»J«»J.'^tOWAL-' 


ORfGfN/\L P,vo:. IS 

OF POOR quality 


- H"*- 


iAC^»^»^A»Ai3LAL)f 

^,^.-feNSI'-^ Svm» i.i(l)» Mi)» Y(l) 

Ilf I i (- ) 

I?0« 1 U i (•< » 

J»SY^l 

Tl ■ I . / ( . / •P'’.4 Cm** 2 ) 
d»C5 l< ) 

CJ 1C l>liJ*UC 
A> AC ( I ) 

G« 1 1 ) 

«1r■2,/t^*A«( r*AAA2^ri**2)/jEN 
Y«»2./x<*j*(T-A**2-C**2)/0tS 
Fh«a^*a*»Y<ay» 

H ■ 0 » 

IP (Wli»'<KNc.O) ^»l./rH 
AZ« -I »•»?(<) -Y*AYZ(K) 
b2»->fT2(K)^vi*x2(K) 

DSl'SOlI + i.- )-bC ( l-lx^ ) 

Al( i) *DS 1 
SZ«C1(‘<)*0S‘^ 

Yl»YC(*^)^ArA.4 2(2.AA*J, f ♦ A* A-f)*0) / AK 

P-b^CJjI PA^YIJ 

C»3YSlMPA*Yi ) 

nGl*.i( i * 1 1 J »L )“G( I~I#JfL) 

naj«Gi i» j^i»L)-5( i» j-i»L) 

OGk»G( I*J»L^ l)“GII»J»L-l) 

L-Al(I)*DafS»*^l(J)*LGJ^CA*XA«>SA*YAA*>/CA(K) 

V'-eKJI^OGJ^SAAAX^e/CACKj-CA^YA 

«»Cl(<)*0&<«^Si*Bl(J)ACGJ^CA*XZ(K)4SA*O^P*Y2U))/CA(K) 

OO-H-^Iaxau-YaavI 

WO»H*( (PAfA^CAAZ)AUAt P*XX^OAFiZ)*V)^OA* 
Nj-MA((P#AZ-JAYX)AUA(PAbZ-JAXX)AV)APAta 
• 0«U0A‘ U V )*‘VO*mO*mO 
SVi n»bIGNIbC'<T(CC)»U) 

IF (ivii»o.6:.o) sv n )-sv( i-i )ASvi i-i )-svi 1 - 2 ) 

OQ»l.AJlA(l,-JOA2.AOMfGAAMA(YXAu+XXAV)) 

Gr( I ) -FMACHASV I I )/Sw«TOd) 

CP( I )-TlA( JQAA3.5-1. ) 
aI»1.A.PASCALA( AAA-JAQ) 

r i»scalaaag 

X(I)»xC(K)>ALOG(SO«T(xIAA2aYI»*2))/Xk 
Y( I )- rC(K) aATAN2{ YI»XI)/X< 

A2( I)»xZ(K}-yxasz 

A3 { I ) • YZ (K )♦ XX*SZ 

10 continue 

kETLPN 

END 


ORIGINAL r 
OF POOR QLmi. I 


-•J J - 


itj*- ji'i iNt c*>Lwr I II. ic*« » jiC»»'“ach) 

L PL.TTS :»» iT tfO*4L lMk-V4L> IN r<»e 

LUrliOcI# »(!)» ‘J C i ) . 

1 « ^ I r ■ c 
.‘»iTc 
po iw 

K LiNkm -KjJt ( i ) 

>-JEN«i./ 0 l U ) 

A A K ■ ) . 
c M A I • , 

Ci>na 

to i-ll. w 

A1A««A''AM(Ai4i.Afj>(l(in) 

DC J J • 1 1 . 

XCBAC-f jes*-*! I ) 

Kl'ItJV./A.iAlMAtJSIKII))**.!, 

L IN£ ( < i ) «K Ji)t ( d I 
K ’ ■ <1 i 

iOJ ) 

LiNc JOc i 3 ) 

1^ vSO TJ 30 

f 3- 1 . 

Gu T O <.J 

30 L iSf- |K j ) •< jnt 12) 

4.0 K<>« 1 

LINE (<<•) ■■C JLt ( 3 ) 

J J ■ 0 

W-ITE I1-MIT.70) XFi AC» * ( I ) » JJfLlNc 
LlNE('^i)»*v30t( n 
L iNfc (K^).K jOE( 1) 

LINE (X3)'K00E( 1) 

LINE ( ] -K DDE I 1 > 

50 CJNTlNlfc 
WETtKN 


4 ^ 1 i\ f 
Jt i ) 


c 


to FQBHAT Cx» /» 3 X» 3 HX/C.^*» 3 H«C 0 ^'i>.<»<» 5 H' 10 SGN. 2 h JN) 
70 FORMAT IIX,F 5 .?.F 9 . j. I 3 , 2 <»iG 0 Ai) 

ENJ 


ORIGINAL PAC:: IS 
OF POOR QUALITY 


i'je->i3jT I Nc (M 

^Otcu AN') THt O'^ofe’fc Cuf^OuFfC Al f ft C -• G - 1 Ti 

PCIM » jHts THE •'aCh *JU**aE*iS A»E w«lTTeN C x IMF C« l^'UI 

G(l?^»^e»l’l.jL(W'i»i7|»c0lin#W(i2;i»i7I.ITtiU7),Iu:(:7) 

l»O(W't)»ftl(12<i)>A2(W'^)#ft3(l?%)*c>0(26)fil(2t>)»'<2(2^)f^3(2r)>7(17) 

2»Cl(i7),C2ll7)»C3(W)»C«»(l7)»C^(l7)»)t:(i7),W(i7),i^2(i7),y;(^7,,y 

2Z(l7l.tZi(l7)»KbY»'#N<.Nr.NZ»'‘TEl»»<Tt2»IS<H»S:AL»S:AL2»*<»J-"fetft»AL*' 

XiHAtCAfSAfFIACHtTf^AfA^LAOE 
OMENiijiN IS (2 A) 
tl» 1./*- 'iAC-«**2*3,2 
A ■ C i ( •< » 

11-2 

I2-N* 

K Y ■ .» Y f i 

0® I T£ ( i J I 

• wilt ( 6» <.C) Mh 

JJ 2 j I*ll»i2 

CSI-^:( ixl.w l-S2( I-i»M 

CS<-S0(4»kM I-SCIIjK-II 

S«» All i ) ftOS I 

SZ«Ci('< l♦0S< 

D3 10 f f<y 

X > At( I ) 

Y-SCI I.-<)^30(J I 

CFs«(Tx<*x-Y*Y )**?^(2,*xftr|A*2 
xx»2./ »i4*x*( T^xftx^YAY )/0EN 
Yt'Z./XWftY^CT-X^X-VftTj/OEN 
FH«x***^yY<a*2 
M»0 . 

IE ( F M • j T • • 1 r - C •? ) H»i,/FH 
AZ«-XXYXZ( K » -Y YZ (K ) 

6Z"-*<ftYZ(<) X*XZ(K) 

Yl»YC(K)^ATAN2(2.AX*Y»Tx-XftX-Y*Y)/X*< 

P»A*C ISiPAftYl) 

0« A*S 1N( FA «Y1 ) 

DGI"G(1y1»J»w.I“&( I“1»J»'<) 

DGJ*0( I»JY l#^)"G( I»J“1»K) 

0Gk*G( 1»J*^yI)“G(I»J»K“1) 

GX«A1 ( I )*DG1 
GY.-bk j »*0«,J 

U»GX-SX»GYYCAftXX4>A*YX*P/CA(K) 

V»GY^SA*>X*P/CA(K)-CA*YX 

w«ci(K)*oG'<-sz*GYx.cAftxz<o^iA*(CYp*YZ(Kn/:A(K) 

00«H*( XX*U-Y X*v ) 

VC«H*((PAYXyO*aZ)*U^(FAXXyQ*3Z)*V)yG*w 
WO«M*( (P*AZ-3*YX)*U4(P*BZ-a*XX 
C«UO*U&AVO*VO^ WO*lxO 


ORiCiN/ L : . 

OF PO i\ » = J w.i 


V 


-93 


J 


1 S ( J ) - SyPT (a / (01-0.2* ( 3-2. •0*ieGA*H* <vx*u**»*vn)l*iccc. 
10 CQNTlNce 

wPlTt (O>:>0) ( 1S( J)*3«2>KY) 

20 continue 

PETUPN 

30 f OP^'AT ( iMi » 

^0 FOPKAT (IBH printout of SP£E0#/»13m FNACh • 

50 FOPNAT (2AI5 ) 

END 


bUbP jul INF F0»CF (II, I2» *» T,CP» al,Cmu«D,»'*,CL ,C0,C^ • *»<,Vk,Z1,Pa) 

CALCULATES SeCTION FQkCE COEFFICIENTS 

DIMENSION <(I), V(l), XK(1), YK(1), CP(il 

WA0«57.2957e 

ALPHA>AL/P AD 

CL-0. 

CO-0. 

CN-C. 

N-I2-1 

00 1C !>11,N 

D>-( M l^l)-x(I ) )/CHOPO 

JY- (Y( lfl)-Y( I ) ) /CM080/ Z I 

XA- ( . 3* ( X ( !♦! I ♦x ( I) )-XH I /CHOP 0 

YA« ,5*( Y( !♦! )^Y( I ) ) 

CPA».5*(CP( !♦! )^C? ( I) ) 

0X-0X*CJS(PA*YAJ4XK(I)*SIN(PA*YA)*0Y 

OY«OY*(l.^YK(n*SIN(PAPYA)**2/2I) 

YA«((1,-PA)»YA^SIN(PA*YA)/ZI) /CHOPO 

DCL --CPAPOX 

CCD»CPA*OY 

CL-CL^DCL 

CO-CO^OCO 

10 C1-CN^0C0*YA-0CLPXA 

DC L -CL • COS ( ALPHA) -CO*S IN (ALPHA) 

C0-CL*SIN( ALPHA )4C0*CaS( ALPHA ) 

CL-OCL 

PETUPN 

END 


t 



UHJGINAL PAGE 15 
')F POOR QUAUTY 


S'». luriNE (KTE » 'C»C3»tL 

If-* , ( » Pa, ArtL 4 CE ) 

C' :o.Afci TOfAL FOk-Ct COEFFICIENTS 

c;' ;i:n :^o^ni)» scud* s:ju)» sc^u)» zu). c^m 

SPA- Z I T t ? ) - Z ( •< T F I ) 

CL-C. 

cn-c. 

C“p-C. 

C l-» • J. 

CFr *0 . 

s»o . 

N»^T£Z”l 

no ij i'»ktei,n 

AZ-. 5 *U./: 5 (^^ 1 MI./C 5 (<) ) 

OZ- 02* ( l.-p A » • ( Z (•^♦i ) -Z ( >< ) ) / 2 . 

AZ»AZ^(l.-PA)* (Z(>^^1)^Z(<) )/2, 

CL»CL»JZ*(jCL('<^l)*CHC<iO(<^i.)^SCL('< » •CMQ.i 0 ( < ) ) 
C0«C0*JZ*(S:0(K^U*CHC-^0(K^i)^SCD(< )*CHUPD( ■< . ) 

cMP-ciiP»oz*(C“«apc(<*-i)*(SCN(K 4 i)*CHO»<o('^>i)-scL('<^n*><c(K*i))*CMC«- 

lD(K)»(SCf('^>*CHOPO(K)-SCL(<)*<C('<))) 

C«R.C^F»aZ*DZ*(SCL(k^1)*CmCrO('<^1)4SCL('^)*ChUFD('<)) 
c'‘Y»c*iT»AZ*nz*(sco(<^i)*C‘-*CFO(*<^i)^sco(K)*c-tOPO(K)) 
S«SoOZ*(CHJ-U(i‘*l)*CHO*^D('^) ) 

A?L AUF “S 
CL-CL/ S 
CC-CD/S 

CMP AS / s **2 

C***' • ( C iK ) / (SAjPAN) 

CMY. ( crrtc ^Y ; / (i*SPAN ) 

P6T UPN 
END 


SUoPJUTiNfc PF.FIN 
HALVES MESH SIZE 

COMMON (,(12=»»26,17),SC(l29,i7),60(17>,lV(12<?,17),nEl(17),ITt2(17) 
l»A0(129),Al(l2'3),A2(12‘/)#A3(129)#u0(26)»')l(26)»«2(£e))»d3(26)»ZU7) 
2#Cl(l7),C2(17),C3(i7)#CA(17E,C5(l7),XC(i7),xZ(17)#>ZZtl7),YC(i7),T 
3Z{17),YZ2(17),KSYM,NX,NY,NZ#'<T£i#*<TE2,ISYM,SCAL,SCALZ»XK,GMtbA,ALP 

aha,ca,sa,fmaCh,t,pa,aelace 

M*»NX^ 1 
KY-NY^-1 
MY-NY >2 
MZ*NZ ^3 
MXO«N* / £♦! 


OP' •. 

ur- eoor< 


fZC-NZ 

L'C 6C 
J«NT/Z^1 
J J • R T 
1C 

1 

2 0 G(Il*Jj**^)*^(l*j>**) 

I- l-l 

II- Il-i' 

IP ( I .or .j) GG TC 20 

J« J-1 
J J • j J-2 

IF (J.oT.O) Gj to 10 
DC JC J»i» X y , ? 

D1 30 I»2»N*,Z 

30 G(Itj*K)>.?«(G(I«l»J»XMo(I-l»j»X)) 

00 5J I-l.'ix 
DC <*0 

<>C G(IfJ>x)>.^»(3(I»J^l*x)^G(l»J-l*0) 

50 G( I fir ,K ) • j. 
cO CONlINUc 
H Z .'< ■ f'. Z 0 
f*ZST-1Z 
70 continue 

00 ao j-i» iY 
OC 90 I-1»'1>1 

ec G( 1#J»«ZST )-G( I »J»MZ^) 

IF (r'ZSi.cO.i) gC to 100 

MZST».XZiT-l 

on yo j«j»'^Y 

00 <iC I-ly'lx 

90 G(I»J»^ZiT)«0.5•(G(!»J»'1Z•*)♦G(I»J»^Z^'-l)) 
H 7 M- -Zf'-l 

IF(MZST.ka.l) GC TO iCO 
f^ZST-'^ZST-l 
GO TO 70 
100 CONTINUE 
KK« 3 

DO 150 X«K<»«Z 
l«f-X0+l 

IF (K.lT.xTEI.CR.k.GT.KTEZ) 30 TO 12C 
Il-ITEI(X) 

I?-ITE2(<) 

DC 110 I«I1»I? 

DSI »S0( I+1»K )-S0( I-1»K ) 

D$X»SO( i #<♦! )-S0( I »K-1 ) 

SY-Ai { n ♦QSi 
SZ-CK K) *OSK 

0GI«G(I^1#KY»k)-G(I-1#KY»K) 

0GK-G( I»KY.K*l )-G( I»KY»K-1) 
R-Af^lNO(l*IV(IfK ) ) 



-Jo- ORIGINAL p;.:,:: t'i 

OF POOR QUALK"^ 

< ■ &C ( I ) 

1 ) 

•• •♦Y**T* 

h • J / A 

aZ«-*»»*« (<)-y<*v/(-< ) 

A « I • C 7 ( ^ ) 

Yl»YC(^MArfi\/(;,*»*Y»T»i*A-Y*Y)/krf 
»-«A*CLi( kA*Yi) 

s« A* WS( kA jri ) 

F //•*.* M 

Fyr«-**^#(Y<«Yi 4 f 2 Z*(r/*'‘Z^»A**X))»l,-- 

(-X«>» 4 «M«(Xt««(«K^Z*(AZ«AZ^YA«YX)) 4 i,-i> 

(AZ*‘^Z^AX*Y< ) ) 

M Z * -t • Z Z • A z 
f y / •-i* 1 1 i * •* i 
i»»F<Y-A»*^»»- 5 Z^'-<A 
AY- FY Y-j A ♦■= » Y-SZYPYZ 
AZ*f-Y?-j>*PY/-SZ*FZZ 

OY-AY-jX’AX-jZYA/ 

J«.M(II» 0 C. I»:A***^SiYY**K/C<-(<>-J''f' 5 AYYA/C<*l'<) 
-«Cl(-')*L,;.^*CA*''Z('X)tbAY(j4>i*YZ{'<))/C‘»('<)-.'t«A*YZ(<)/C^l‘<) 

lie Cl I.KT-i»*«)^(AA*jYAY*V^AZY-) / (tiY^jllKY)) 

EC(’a)*'j( lt*^Y»K ) — ij(ll»^Y»(^ ) 

1 »1 1 

li?0 1«I-1 

c *0 t 

IF ( I j ( 1 »< ) . . 1 ) GJ T J 133 

F »i- W ( A ) 

13C 

(jll»AY^l,K)«.i(*A»KY-i.K)-C 

GI'1^<Y♦1,^).^( I»KY-l»K)tt 
IF (iVlif'xj.Nc*”!) GG TJ l^^j 

Gll»Ay,<).,‘>* 3 (I,KY»K-l)^. 23 *(&ll»KY»K^i)^G('^»AY,-^ 4 l)) 

IF IIv(i»A«.l).LT.i) G(I»aY»M»,!j*G(I»ky, <♦!)♦, 25♦(G(I»<Y»^-i)♦G|A. 

1KY»K-1 ) ) 

G(M»^Y,l<)«GII»»<Y,-() 

GdiKY-l^Kj-.^^IGII.KYjO+Gll.KY-Zi-O) 

GIM,^Y-l»0«.3*{G('^fKY.<)4'G( "AjK Y-2» K ) ) 

1<»C IF I I ,oT .2 ) G J TC 12C 
15C CCNTlSv^t 
f<ETU» >1 
E^0 



cr?iGif. 

-Q7- OF POu.C I,. 


C cl'oiC a*'LlNP »=IT •'*^SC**i^E0 CV«JIT1C'! 

LIflfcoiJN j(i)i Mi), •'.^(l), FPJ(i), hi>->a(;) 

( ric.ko.o) &c rc !»• j 

1 s u ■ C 

IF ('^“1) l^'».l^C»iO 
1C K- ( r -»^) 

I-'i 

J ■ ♦R 

O ^ ■ S ) J ) - i ( I ) 

O-Oi 

IF (OS) ?t.irtj.20 
^C UF • ( F ( J )-r ( I ) ) / OS 
IF ("i-i;) jC.^C.SC 
3C U-.5 

V*i.*(t»F-v'1)/0S 
GO TO 30 
<• C U ■ J . 

V ■ V ^ 

G t T w 0 0 
iC U»-i. 

V«-GS* V 1 

cn n co 

60 I ■ J 

J ■ J ♦!< 

C S • S 1 J ) - b ( I ) 

IF U*bj) i?C»l?C»70 
70 CF- (F ( j)-F ( I ) ) /CS 
3»1./(LS+0S^U) 

U«-J*OS 

V»i* ( b.»CF-'/ ) 

80 F P( I ) -0 
FPP { I ) « V 
U ■ ( I . -0 ) *0 S 
V»6 , *0F^ DS *V 
IF (J-is) o)»^J»6C 
<?0 IF (^N-£) ILC.11C.120 

ICC V ■ ( fc . » ) / U 

GO TO liC 

lie V-VN 

g: tc lie 

12C V-(0S*vN 4F-»P(1 ) )/ (l.*F-»(l) ) 

130 e-v 

o-os 

1<*0 CS-S< J )-S>( I ) 

U»FPP(])-FP(I)*v 


- )«1- 


ORIGINAL PAGt IS 
OF POOR QUALITY 


PPP( I ) *0 

V-0 
J« 1 

I .I-K 

IF 1<»0» I iO» 1<*0 

liC l-N-" 

FPP?(>0«FP‘>P(I) 

FPP (N ) • d 

FP(N)«0*-4L/*(FPPII)^3^t')/''. 

1N0« 1 
IF (VJLt) 

IfcC FPFP{J)»ri-1 
V«FFP ( J ) 

17C I-J 
J » J 

U«F PP ( J ) 

FPPP(J)■^•PPP(I)♦.3♦C-S♦(F(l)♦F(J)-0S♦^i•lJ♦V)/l^.) 

v»u 

IF (J-M 17J»1«C»17G 
18C Continue 
F C T L X N 
END 


SU3><DoTiNt INTPL ('1I»NI»SI»FI»h»N»S#F»F»'»F?P#FP?P»^3LE) 
1NT6*«PCLAT IJN OF CLsIC iPLiNF 3Y T&YLC? SEPlti 
DI'*ENiIGN >I(1)» FI(l)» S(D# F(l), FP(i)» -PPd)# rPPP(i) 
K« I aBS ( N-'* » 

K« IN-M )/K 

I N « 1 
N I » N I 
0«i I ) 

IF I0»(3l (NI )-SI(^I ) ) ) 10»ZG»2C 
10 ^IN«NI 
N IN »f1I 

20 KI- lAdb(MN-r*IM 
IF (KI) <»0,A0»3 0 
30 Kl- (MN-MIN) /K I 
40 II-KIN-KI 

C-0. 

IF (NODE) oO»30»50 
bC C-1. 


SS-SIllI ) 

70 

(l”N) 

90 IP OMi(I)-iS)) 7C»7('.VJ 
9C J-I 
I-I-i' 

SS- SS-i( I ) 

FF«FPP^( l)*.23*Si^FFPPP 
FF-FPP (1 )♦ iS*PF/3. 

FF -FP ( * ) ♦ . b j»FF 
FI(in»*^ll)^ji*FP 
IF (II-NlN) 

ICC Cr-NTlNLc 
•^ETU-’N 

e^ D 


SUBFoUII'<fe 1 "<^Ecli (IPLJT»Sv»ii*1»CP|*»7»rirLE»OC»AL»2JNt»P''''2»jEv»CH3 
IkJ0> XiC'tL>PSCAL ) 

GFNt-rtTcj THSrE C I HE N 5> I'JN A l C&LCOHP PlTTj C^C 9c0o 

G(W-^.>b»17)»SC(lP7»i7)#E0(l7)»IV(12^»17),n£l(in.ITE2n7) 
l#A0(12t|,A4(i2<J),A2(12^),A3ll29)»3C(26)»U(2C5)»r?{26).93(2fc)»2(I7) 
2»Cl(i7)»C2(l7)»C3(17),C^(i7),C5(l7).xC(i7),Ji2ll7),XZ2(l7).yC(i7),Y 
32(17)»TZ2(l7).K5rH,Nx,Mf,NZ#KTti»<TE2»IjfH,5CAL.SCALZ»x<»0''tGA,ALP 
^tMAfCA* i»A#F iACi » T #PA» AfLAJE 

OlicNilJN x(l), vm, S\/U), SH(1), CF(I), riFLEdC)* P(2 j) 

M«1 

IP ( XSC AL ..-ic .0 . ) SCALX..3 *AbS(XSCAL)/ChJPOC»ZOnE 
DZ»Z(KTc2)-Z(<TEl) 

IF (PA.ci.l.) C■Z»l./C5(^rF2)-l./C5(<TE 1) 

IF ( PSC AL. 3C . J . ) SCALX»^./0Z 
SCALP«-i .0 J 

IF (PaC AL. 4E .0. ) scalp--. 5/AdS(PSCAL) 

T X » 3 . J 

SX»-SCAl>*xC (MEl ) 

IF (IPLGT. NE.il GO TJ 10 

CALL PLOTSdL ( 1 OuO » 2 5 H AN T Q I NE duURGEAOF ♦337»^H) 

1C continue 
IPLOT-C 
call FPAHF 

CALL plot (1.2!?»I.»-3) 

ENCODE (ij!?»l<^C»#P) FmaCM#F;12,0EV»AL 
call SYHbOL (C.C#C.73# .lA»Rf G.,65) 

ENCODE (oC»200aR) 


20 

30 


4.0 


5C 


-1 o:- 


COMIM »: 
C?NT iNwt 



K.i 

C JM I I NL'f: 

■ n ♦ i 

u (•< .or .< rf 2 ) GO ro 7w 
IF (K.Ll.^rtl) GC Tj :J 
Il-ITEi(K) 

I2■Irc^(^) 

CAlL </tLL 

bC3«3.Mi(<)-Z(<TFl))/C(<rt2)-2(*«Ttl))^2.7S 
DO ‘♦C I» 11 » U 
MI 

Y(I)«5LALX*MI)tSY 
CP( I)-SCal?*CP( IMSCP 
CON 1 I '.Lc 

IF (r.Lw.i) Cj to t-0 
N»I2-I !♦! 

call LiNt (<(ll)»CP(Il)*Ntl»C#2>o.fl*»J*»i.) 

GO TG iC 

continue 


N«I2-Ii^* 

OC ^0 1« II » 1 2 
x(i)-*(i)4>rx 
bC CGn r i xOc 

CALL LiN( (X(Il)«f(Ii)»N*liJ»2fJ.fI.*C**I.) 
GC TJ 2C 
70 CCNTiNUt 


IF (M.GT.i) GO TO 80 
GO TO cJ 
80 CONTINUE 
CALL Fi'Art 

CALL ^tLC (*<TEl»KTtl»SV»S'l»CP»x#Y) 
Il-ITEi(KTEl) 

I?» ITE2(kTE1 ) 

SCAL x»l , ii»XK /3 . i<.l8926a 


UG 90 I-I1M2 
X ( 1 ) «2 . *x ( I ) ♦SC AL X 4SX + 1 , 
Y(I)«2 .*T(I)*SCalx 4,5 
CP ( n» Y ( I)4b. 

90 CCNTINL6 
N-I2-I !♦! 

CALL LINE ( X ( I I ) » Y ( I I ) > i *Uf 2*0. > 1 .f C . « I . ) 
CALL LInE (x<ll)fCP(Il)>N*l*0*2*0.*I.f0.*I.) 
ENCODE (6O*210*R) 

CALL SYMBOL ( I .C»3.?» . 1A,P ,0. » tO) 

ENCODE (bO»22J»P) OC 

CALL SY«bOL (1.0*7. 75*. iA,R,U.,bO) 

CALL PLOT (-1.25*-l.*-3) 


ORIGM'iA.t. i t t J 
OF POOR QUALliV 


R y ■ Nr ♦ i 

N A ■ P A 

K • « T t i 

ICC CJN1 : 4UC 

call r-^AM 
call •>LuT ( 3 . 

0<3 lU J«C»<Y 

u«3L(J)^bw(L<iK) 

A■*C(^)♦AL3J(A•jS(l.-.l.•iv.AL»^♦ll))/;< 
IP ( J*'.j.L T. r ) G.T TU n: 

110 C"?nTI''<Lc 
ICC CONTiNLt 

o-eo(N T ) ♦s ) ( c . •« ) 

9-l.^U-(C)*Aj(C )-l*2 ) / r 
J»SCal*al(< 2 )*j 

A »-2 . / A 
H ■ » / 3 

IF (A.tc.C*) A «3 
SSa» 1 .C / 

SS<»Aii\l(%S*»A» 6 ) 

DO 130 i«r»N< 

A« AO ( I ) 

L®L 0 T -3 
DO 13 0 J «?*•<» 

C»oC(J)^aO( l»N) 

Xl -1 .♦^•SCAL* ( A*A-C*L) 

Vl»iCAL*A* J 

X 2 »*C('^)+ALL‘Olxl*Al 4 ri*rL)/AK/ 2 , 
Y 2 »YC(K )^ATAV. 2 (Y 1 »X 1 )/XK 
Y 3 » 3 IN(Y 2 )/C 3 (K) 

X 2 »SiX*X 2 
Y 3 »SiX*T j 

call plot ( X i ,Y 3 »LPLJT) 

L PLCT -2 
130 COnTHUc 
LPLOT- i 
DO 1<*0 J -2 »KY 
C« 30 ( J)>iO(LX, K ) 

Xl- 1 .-. 3 *SCAL*C*G 
x 2 »xC(<)+AlOG(A 3 S(X 1 ))/xn 
Y 2 « YC ( K )-A TANf ( C. . X 1 ) / XK 
Y 3 -SIN( Y 2 )/C 5 (K ) 

X?«bbx»x 2 
Y 3 » iSx *Y 3 

CALL PLOT (X 2 » Y 3 »LPL 0 T) 

LPLOT -2 
1 <.C CONllNofc 

DC IbO J« 2 »KY 


ORIGINAL r 
OF POOR Qj/ 


I plct. j 

C!l lou !■£»'<* 

A*AC( 1 ) 

A* 4-J*0 ) 

r 1«5 )Cac»a* 1 

V2» y: (m )^a ta\> ( yi , »i 1/ »< 

Y3»S1N(Y^)/Cj(*^l 
X?- « • < t 
Y3"Si**Y3 

IF ( l.^.L .L« ) jC 10 15C 
Y2 -?.*tC I*' l-Yi 

Y<.»S:)**iI.NlY2)/Cii(K) 

call PL j I U J » Y^ # LPL jT ) 

LOLoT- 3 
li»0 CInTINoc 

call PlJI 1 » ?» Yi, lPL jT ) 

LPLGT-2 
16C CTMlNOt 

i 2 n ^ M*n<) 

ENCGJc 

CAlI SY'^tGL ("l»C»A.t»«lA»i^>C»»^J) 

^all PljT ( “ 3 » » -a , 5 » - 3 ) 

•« 0 

IP K.Lfc.'<Tf2) GO Tj lOJ 

K aK-l 

IF K.Fw.<r£2) GC TO 100 
CC» 1 . / YK 

IF (PA.Nt.l.O) GO TO leo 
N«2 

CO 17u 1-i.N 
S-FL0 aT( I ) /F L JAT( \) 

X0«0.5/ZuNc*S 

call C'jT («0»i')C»SCALP#jY»S«»CP»A»Y»Z»TC.ITEl»lTt2»*<TEl,<Tfc2»-'SY'<) 
17C CONTlNUt 
IbC CjNTlMCt 
fcETUSN 


19C 

FCP EAT 


16HALP ■ 

200 

format 

210 

fcpmat 

220 

FC'PMAT 

230 

FOPE AT 


END 


(5H11 ■ »F<*,2*1h» 2*»3HM2 ■ »Ff«2#lH»2»#0^nEV 
»FA. 1 I 

(ildPPESSjPE 0IiT«<IeuT If>^»5i( » Ui3L AOE PrfjFiLt 
(t3M CASCADE PEPPESENTAT ICM 
(9H G/C »»F5.2) 

(25H GPID ON The jUPEACE i •»Fj.2) 




-ijj- OF P( ‘.* 


l^ ‘ r ^ ) 

C Tr-Ii Swj^OjTl^E *»l:TS SiCTIjN3 JF Tnr C ‘j i ► ^ i!.c 

CC'I^JS /iw**/ «1( I ti * 2 t l*G^(l2^»26)»S*(le^)*^t •/( i,' 

t)I^*t•oiJ^ ivu), ^'i(l). :p(i)» Ml)» y(l), /(i). iTt.'Ci# 

iYC ( 1 ) 

Call 

NZ-'^r' 

CC so •‘•'•Tri.^TEZ 
call vtLu 
i 1 o ■ 1 T ‘ L ( •< ) 

I ,.'C • I r t ) 

e« 1 • 

CO n i-iii»i2c 
IF ( < ( I ) .t r , ic ) G J T j £‘J 
10 COM INLC 
I- I-l 

R-C . 

?C COM I ilE 

F • ( *u-M 1 ) ) / ( M I-l )-A( n ) 

SV(«)»-*(r(I)^F*(Y(l-i»-y(i»)-YC('<))4YC(<) 

C(")«(t''auFMCi»(i-i)-Ci’(nj)*^*scAL»'-z. 

s ■ I . 

on i J !• ii )» Ww 
M- 1 1-> ii:-i 

IF ( < ( '1) .L r . X J ) oJ n ^0 

30 COnTI^lF 

o«C . 

<»0 CONTI yL't 

F •( XO-X ( ," ))/{*(-♦! )-x ( .1 ) ) 

'J(^)•{C»»(^^♦F•(CP(f^♦l)-C‘>(^1)))♦3•SCALF♦•^. 

p^(^)«EAP(Z(><)-Z(<Tk2)) 

«!(•<)«?. ♦*<0(K)-l, 3 

50 CONFlNUc 

CALL PLOT { l.f 5.»-3) 

Call LIU (•'i(KTkl)»C(^TEl)fNZ*lfO>£io.>l.»0.»l.) 

CALL LINE (•<l(MEl)pl3(^rcl)»NZ»i»0»2»C.»l. 

K*l./DC 

CALL PLOT (X). » C . »-3) 

O'J fcU L*l#'^ 

XL-L 

DC to K-KrSl»-^TE2 

X(K).j.A*.:(O*C05(SV(*<)^2.*A.lAi5v265AXu*DC) 
Y(K)-3.*FO(K)*SlN(Sv(K)4-?.*j.l<.l59 2b5»XL*OC) 

to continue 

call line (x(<Tci)*Y(KTEl)»NZ>l»C»2*0.>i.f ). »1.) 

00 70 K-KTE1»KTF2 

A«2.a J,i<»l392‘55*XL*0C 

X(K).i,A»,0(K)*COS(Sf*(K)*A) 


ORIGINAL PAG2 ’V’ 
OF POOR QLALHV 


Y(« s^l" I^A) 

7C COf^UNUc 

CALL Ll'<t (>i(<TYi)*Y(HT£l)fN^#l>Oft!*C«(4*fO*»lt) 
^ C C 1 i «ot 

A ■ 3 . • ' « »* ( 2 ( " I ; i ) - n X Tf ^ » J 
03 ^C* I"i»Y7 

T'i.^i.lAU'Ji^tAFLJAn’-U/V^. 

Cl I l•^.•C Jil I ) 

Cl I » «3 . *5 1 »( n 
X I 1 I • A*CuSI T ) 

Vll»«A*jiMM 
90 CONTlM't 

Call Ll^ir |X{l)tYli)»97fi»G»2fC**It*\/*fl(l 

If l-'ST''.?;.!.) C3 rj ioj 

call Li^ib icii)>cm>97»i#o#2»r.*i.to.*i.i 
IOC CQNTInlc 

XJa3.*Xv^/ A 

ENCOOt loJ»U3»*<l) x: 

call SY.idJL i-2.i*A.;*.:A»KA*c.*tc) 

CALL »*LJl 1-7. #-5. #-3) 

FeTo*^ < 


lie f3-<iAr livi SECTION IN TrtE PLANc x •»F5.3) 
ENO 


